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Executive Summary 


A Delta II rocket exploded seconds after liftoff from Cape Canaveral Air Force Station (CC^S) on 17 
January 1997 The cloud produced by the explosion provided an opportunity to evaluate the models which are 
used to track potentially toxic dispersing plumes and clouds at CCAFS. The pn^ goal of this Project was to 
conduct a case study of the dispersing cloud and the models used to predict the dispersion resulting fi-om ± 
explosion. The case study was conducted by comparing mesoscale and dispersion model results with available 
meteorological and plume observations. This study was funded by KSC under Applied Meteorology Unit 
(AMU) option hours. 

The models used in the study are part of the Eastern Range Dispersion Assessment System (ERDAS) and 
include the Regional Atmospheric Modeling System (RAMS), HYbrid Particle And Concentratmn Tr^sport 
(ElYPACT), and Rocket Exhaust Effluent Dispersion Model (REEDM). Two different versions of RAMS were 
used in this study: RAMS version 3a from the ERDAS prototype system and RAMS version 4a from the 
Parallelized RAMS Operational Weather Simulation System (PROWESS). 

The primary observations used for explosion cloud verification of the study were from the National 
Weather Service’s Weather Surveillance Radar 1988-Doppler (WSR-88D). Radar reflectivity measurements of 
the resulting cloud provided good estimates of the location and dimensions of the cloud over a four-hour penod 

after the explosion. 

Meteorological data from local observations and sensors provided a basis for comparison with 
meteorological model output. The meteorological condiUons on this day were strongly Muenced by synoptic 
rather than local forcing. Observed data were obtained from the WINDS tower network at CC^S/Kennedy 
Space Center (KSC), rawinsonde data from CCAFS, 915 MHz and 50 MHz radar wmd profilers at 
CCAFS/KSC, and standard local station and buoy observations. 


In this report the terms plume and cloud are used interchangeably. In dispersion meteorology, a plume is 
usually thought of as an elongated mass of gaseous or particulate material that emanates from a source and 
widens and disperses with the flow of the wind. A cloud is usually referred to as visible mass of gas or particles 
sh^d like a sphere. When the Delta II exploded, it initially produced a large cloud that dispersed and moved 
downwind in a plume-like fashion. 

The conclusions of this smdy can be categorized according to the plume observation technique and 
according to the models used in the analysis. The findings of this study are: 


WSR-88D radar as a plume observation tool 

• The WSR-88D is a good tool for providing plume t r acks from rocket explosion plumgg . The radar 
provided excellent data on the location and track of the resulting potentially toxic cloud. The data 
were extremely useful for model verification since no ongoing program is m place to me^ure 
plume track or concentrations. Bud Parks of ACTA, Inc. is conducting a study to capture data from 
nominal and abort launch clouds and has captured radar data from three Eastern Range abort 
clouds and 36 of 46 (78%) nominal launch clouds (Parks and Evans 1998). 

• The WSR-88D does not provide concentration data . The only data obtained by the radar is die 
reflectivity value measured in dBZ. While this measurement gives an estimate on the relative 
density of material (smoke particles, water, dust, or other particulate matter), a methodology is 
needed to convert dBZ to concentrations of hydrochloric acid (HCl), nitrogen tetroxide (N2O4), or 
other materials of interest. One of Range Safety’s main concerns is determining the exposure limit 
(concentration over a specified time) of certain toxic materials. 

A dark orange cloud at the very top of the large lower cloud was initially visible. The dark orange 
cloud most likely contained some amount of nitrogen tetroxide. B^ause it was located near the top 
of the cloud, it is unsure how much, if any, of the N2O4 mixed within the cloud and made it to the 
surface. Our analysis was not able to determine the concentration of N2O4 in the explosion cloud. 



• Vftrtiral plnmpi height data for this case was not very accurate . The radar appeared to accurately 
track the clouds’ trajectories in the x-y dimension. However, the vertical measurements appeared 
to be inaccurate for two reasons. The first reason was that for long distances from the radar, such 
as the 35+ kilometers distance from Melbourne to Cape Canaveral, the radar beam widened 
enough to introduce inaccuracies in the vertical plume height measurements. The second reason for 
inaccurate vertical measurements was because of the strong inversion causing the radar beam to 
bend and bring about measurement inaccuracies. 

RAMS model 

• The vertical wind profile predictions in this case show agreement with observations. Both ERDAS 
and PROWESS configurations of RAMS produced wind flow measurements that matched closely 
with rawinsonde and profiler measurements and seemed to provide adequate input to HYP ACT . 

• RAMS under-predicted onshore flow at the level of the Delta II cloud . Both versions of RAMS 
predicted onshore flow in the 600- to 900-meter layer in the area south of Cape Canaveral. 
However, the observed winds had more of a northeast, onshore component than the RAMS- 
predicted winds. The movement of the actual explosion cloud tracked onshore north of the model- 
predicted cloud. The prototype ERDAS RAMS predicted winds with more of an onshore 
component than PROWESS RAMS. 

• RAMS accurately predicted the strength and height of inversion for this case. The well-defined 
inversion that was measured by rawinsonde and based at approximately 800 meters above the 
ground was predicted by RAMS to be based at approximately 750 meters above the ground. The 
inversion, as determined by the vertical temperature profile, had a significant influence on the 
explosion cloud. 

REEDM source term 

• Characterizing the source term of umque explosions is difficult. If a rocket explosion occurs, the 
circumstances will be different each time it happens. For example: What was the flight time? How 
much fuel was consumed? What were the height, location and distribution of the explosion 
products? Were the hazardous and toxic materials separated or mixed within the cloud? Did the 
second stage ascend and then explode as with the Delta II? Did the solid rocket motors explode 
immediately or did they follow an errant path before they exploded? All of these questions make it 
difficult to develop a model that will accurately assess and characterize the source term. We were 
able to use information obtained after-the-fact from radar and video to characterize the source term 
but in a real-time scenario only estimates of the source term characterization can be made. 

• Splitting the source into two sources for HYP ACT model was a reasonable approximation, In the 
post-explosion mode of this analysis, we split the source into two sources for the HYP ACT model 
based on observations. This methodology proved better for this case as opposed to using the single 
column source term that was generated by REEDM, In real-time, REEDM currently does not split 
the source term. 

HYP ACT model 

• Plnmft onshore further to the north and earlier than predicted. H YP ACT moved the large 
lower cloud resulting from the Delta II explosion onshore at a point that was ^proximately 12 
kilometers south of where it actually came onshore. HYP ACT predicted the plume would come 
across the coastline in the Satellite Beach/Indian Harbor Beach area when it actually crossed the 
coastline in the Cocoa Beach area. The observed winds at the level moving the cloud at 600-900 m 
were northerly prior to the explosion but then shifted to northeasterly during the 1-2 hours 
following the explosion. RAMS predicted the shift of these winds from northerly to northeasterly 
but did not predict the shift as quickly as observed. RAMS supplies wind and meteorological data 
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to HYPACT. Because of RAMS’ gradual response to shift the winds, HYPACT missed the 
location and timing of the plume impact on the coastline. 

Tr.i.r.torv diffusion and timing of HYP ACT plumes showed sinularities to observed pl^ 
Except for the problem mentioned above, the trajectory, diffusion, and the timmg of the HYPA 
plumL were similar to the plumes observed by radar. One favorable result was noted m the spread 
and diffusion of the lower cloud as it moved south. The cloud spread m the crosswmd direction a 
a rate and in distance similar to observed. 

Ranne Spfrtv’ s RF.EDM predicte d the movement of plume to the south (176 degreesL .^EDM 
using the 1613 UTC rawinsonde fi-om Cape Canaveral moved the plume to the south and kep i 
offshore until it reached the Melbourne Beach area. REEDM did not account for the wmds with ^ 
easterly component that existed at a height of 700-800 meters in the area over the ocean to the 

south of Cape Canaveral. 


Recommendations 


De velop methodology to correlate conc e ntrations with radar reflectivity measuremente ^The WSR- 
88D proved to be a valuable tool in tracking nominal and abort rocket plumes. However, the radar 
provides no information on the concentrations within the clouds. What is needed is measurement 
of concentrations within the plumes using a sample collection method or another remote sensing 
technique such as lidar. This data could then be correlated with radar measurements of reflectivity 

in dBZ. 

Tmnrnvements are needed in HYPACT plume dynam ics algorithins^HYPACT currently treats 
plumes as non-buoyant, non-depositing entities. We recommend that future enhancements should 
be made to HYPACT to improve its ability to handle buoyant plumes and particle deposition. 
These improvements would allow HYPACT to model rocket exhaust plumes better than the 
current version of HYPACT. 

CnnHuct other studies of rocket explo sion plumgs^Since the explosion of the Delta II two o&er 
rockets have exploded after launch from Cape Canaveral-Titan IV on 12 ^ 

III on 26 August 1998. In both cases the explosion clouds were tracked by WSR-88D racto. 
Detailed studies should be conducted to verify mesoscale models, diffusion models, and radar 
tracking techniques. 
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1. Introducdon 

The primary goal of this study was to conduct a case study of the dispersing plume and cloud resulting 
from the Delta II explosion on 17 January 1997 at Cape Canaveral Air Force Station (CCAFS). The case study 
was conducted by comparing mesoscale and dispersion model results with available meteorological and plume 
observations. This study was funded by KSC under Applied Meteorology Unit (AMU) option hours. 

The models used in the study are part of the Eastern Range Dispersion Assessment System (ERDAS). 
These models include: 

• Regional Atmospheric Modeling System (RAMS) - A nested primitive equation mesoscale 
meteorological prediction model. 

• HYbrid Particle And Concentration Transport (HYPACT) model - A Lagrangian particle 
diffusion model. 

• Rocket Exhaust Effluent Dispersion Model (REEDM). 

A description of ERDAS and its models is found in Lyons and Tremback (1994), Evans (1996), and Evans 
etal.(1996). 

The primary observations used for plume verification of the study were from the National Weather 
Service’s Weather Surveillance Radar 1988-Doppler (WSR-88D). Radar reflectivity measurements of the 
resulting cloud provided good estimates of the location and dimensions of the cloud over a four-hour period 
after the explosion. 

Meteorological data from local observations and sensors provided a basis for comparison with 
meteorological model output. The meteorological conditions on this day were strongly influenced by synoptic 
rather than local forcing. Observed data were obtained from the WINDS tower network at CCAFS/Kemiedy 
Space Center (KSC), rawinsonde data from CCAFS, 915 MHz and 50 MHz radar wind profilers at 
CCAFS/KSC, and standard local station and buoy observations. 

This report describes the modeling analysis which was conducted using the ERDAS models. Section 2 
provides background information on the ERDAS system, the observed radar and meteorological data, the 
explosion and provides some references to analyses conducted by others on the Delta II explosion. Section 3 
presents results of the radar data reduction effort and the mesoscale and dispersion modeling analyses. Section 
4 presents conclusions of the study and recommendations for future operational analyses. 
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2. Background and Study Procedures 


2.1. Delta II Explosion Event 

The Delta II rocket was launched from Launch Complex 17 (LC-17) at CCAFS at 1628 UTC on 17 
January 1997. It exploded 12.5 seconds after liftoff at a hei^t of approximately 438 meters. Figure 1 shows a 
photograph of the explosion. This initial explosion destroyed only the first stage and the boosters and produced 
a large cloud extending from the ground upward. The Delta ff is a three-stage liquid-propellant vehicle with 
nine solid-propellant strap-on booster motors. The second and third stages and payload survived the initial 
explosion and continued upward to about 760 meters at 22.4 seconds. Destruct signals were sent at this point, 
and the exploding second-stage formed a buoyant cloud that rose above and separated from the lower cloud! 
Detailed mformation on the rocket and explosion scenario is found in Ha and Deane (1998). An example of 
damage caused by the falling debris is presented in Figure 2. 



Figure 1. Photograph of Delta II explosion at 1628 UTC on 17 January 1997. 


The clouds drifted in two primary directions. The lower cloud drifted to the south over the next few hours 
passing ovct the Atlantic Ocean south of Cape Canaveral, over Cocoa Beach and Melbourne and continued south 
through Indian River County. The upper cloud moved toward the east remaining over the ocean. The cause of the 
formation of two separate plumes was strong wind shear across a temperature inversion at about 910 m. Winds 
below the inversion were on the order of 1 1 m sec ' from the north and north-northwest. Immediately above the 
base of the inversion the observed winds were within a few degrees of due northerly, then backed rapidly with 
height, and became northwesterly above 1350 m and westerly above about 2400 meters. Descriptions of the 
meteorology and plume track for this event are discussed in Pace et al. (1998), Parks and Evans (1998) and Ha and 
Deane (1998). 
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Figure 2. Damage crater resulting from falling debris from Delta II on 17 January 1997. 


2.2. Data Sources 
2.2.1. Radar data 

The Weather Surveillance Radar 1988-Doppler (WSR-88D), also known as NEXt generation RADar 
(NEXRAD), is a powerful, modem radar wiUi advanced computer processing to provide a wealth of 
atmospheric data with high accuracy and resolution. The rad^ operates at an S-band wavelength of 
approximately 10.7 centimeters with an operating power of 750 kilowatts and a peak power of 1 megawatt. The 
WSR-88D antenna is a parabolic dish 28 feet (8.5 m) in diameter, which produces a 1 -degree diameter conical 
beam. The antenna rotates continuously in azimuth at a maximum speed of 5 rotations per minute and moves in 
predetermined incremental steps in elevation from 0.5 to 19.5 degrees. The elevation steps are determined by 
the volume coverage pattern (VCP) being used. The radar can operate in storm-mode to detect precipitation or 
in clear-air mode when no precipitation is within range. In storm mode the antenna turns rapidly and completes 
a full volume scan cycle in 5 or 6 minutes. In clear-air mode the radar antenna turns more slowly and transmits 
a more powerful pulse of energy to detect targets weaker than precipitation. In this mode the volume scan 
pattern is limited to scans at 0.5, 1.5, 2.5, 3.5, and 4.5-degree elevation angles. The antenna rotates frill circle in 
the horizontal for approximately 2 minutes at each elevation angle before stepping up to scan the next higher 
angle, completing the full volume scan cycle in 10 minutes. 

It has been known for some time that the WSR-88D can detect a wide variety of phenomena. These include 
raindrops, dust, insects, birds, smoke, and refractive gradients. Other studies have used Doppler radars to detect 
smoke plumes. Rogers and Brown (1997) used an S-band (10-cm wavelength) scanning radar to track the 
advection and spread of a smoke plume generated by a large industrial fire in Montreal, Canada. They 
concluded that radar observations may be useful in monitoring plume behavior and estimating the amount of 
particulate material in plumes, although such estimates require better knowledge of the refractive index of die 
material than now seems to be available. Banta et al (1992) observed smoke plumes from forest fires using 
Doppler radar. 


3 


It wasn t known how well the radar could detect the abort cloud from a launch vehicle failure until the 
Delta II explosion at Cape Canaveral was observed on 17 January 1997. A detailed analysis of the Delta II 
abort cloud will be described in Section 3.1.4. 

At the time of the explosion, the National Weather Service in Melbourne monitored the track of the plume 
using the WSR-88D. The Melbourne WSR-88D radar is the closest one to Cape Canaveral. Images obtained 
from their archives are presented in Figures 3 and 4. These images show how clearly the radar detects the 
smoke, water, and other particles contained in the cloud. Local television stations also tracked the movement of 
the cloud using their weather radar and then broadcast these images to the public. 
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Figure 3. WSR-88D image from National Weather Service captured at 1625-1635 UTC on 17 January 1997. 
Launch occurred at 1628 UTC. 
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Figure 4. WSR-88D image from National Weather Service captured at 1655-1706 UTC on 17 January 1997. 


The WSR-88D was a very good tool for tracking the location of the Delta II abort cloud. However, it was 
not able to provide concentrations of the toxic emissions resulting from the abort. The Delta launch complex is 
located 37 km from the radar antenna at an azimuth of 14 degrees from north. On 17 January, the radar was 
operating in clear-air mode and had begun a new volume scan series at 1626 UTC. The lowest elevation scan at 
0.5 degrees would have ended and the 1. 5-degree elevation scan just begun when the Delta II rocket lifted off at 
1628 UTC. Therefore, the 0.5-degree scan showed no return from the LC-17 area. However, the 1.5-degree 
scan detected a strong echo directly over LC-17. The radar beam centerline was 1 123 m above the ground and 
the beam was more than 610 m in diameter at this point. The strength of the echo was 46 dBZ as determined by 
the Level-2 data archive. 

The radar also could not provide a precise height of the base or top of the cloud. An error on the order of 
600 meters can result from the diameter of the beam at these ranges. An echo-free elevation slice above and 
below the cloud is also needed to bound the error. No slice was ever completely under the cloud, and with the 
maximum scan angle only 4.5 degrees, the cloud tops were only estimated as well. An even greater error in 
cloud altitude is likely because of the 5 degree Celsius temperature inversion near 900 meters. The height of the 
radar beam is calculated on the basis of normal atmospheric refraction of the beam. However, when 
temperature increases and humidity decreases as in this inversion, the radar beam is bent downward from if s 
normal path (Doviak and Zmic 1993). Therefore, all radar measurements of cloud heights above 900 m are 
almost certainly too high by a significant amount due to superrefraction of the radar beam. 

After the explosion, the Melbourne WSR-88D data for 17 January were obtained for analysis. The analyses 
were conducted using the software package 88Display (Hoffert and Pearce 1996). 88Display contains features 
which allow users to display detailed Level II data and to obtain various parameters associated with the radar 
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echoes. Level II data are base data which include base reflectivity, mean radial velocity, and spectrum width. 

The features of 88Display include displaying composite reflectivity for all elevation angles, display base 
reflectivity for each elevation angle, building animated time loops, and also displaying vertical cross-sections. 
The cross sections can be built using any two user-selected points from a map of the composite or base 
reflectivities. Another very useful feature of the software is the information query flmction. When a user clicks 
the mouse on an echo feature on the map, the program returns data on elevation, azimuth, ground range 
(distance from radar), latitude, longitude, and height above ground. This information is based on elevation 
angle of the currently selected scan. 

The data collected by the WSR-74C located at Patrick AFB were not available for these analyses (Herring 
1999). The tape containing the data was produced using the McGill processor and its associated software. 
When the WSR-74C was upgraded to the Interactive Radar Information System (IRIS) SIGMET software in 
1 997, the new system would not read the tapes created on the older machine. Having no requirement to retain 
these tapes, the contractor recycled them for other use. 

2.2.2. Meteorological data 

The meteorological data used for this study consisted of the data which are routinely available within the 
Meteorological Interactive Data Display System (MIDDS). The RAMS model runs which were made for 1200 
UTC on 17 January 1997 were initialized with the following data: 

• 1200 UTC Weather Information Network Display System (WINDS) tower data (39 stations) 

• 1200 UTC surface and buoy observations (484 stations) 

• 1200 UTC rawinsonde data (37 stations) 

• NCEP Eta grids generated from the analysis at 1200 UTC 17 January 1998 for the following 
forecast times: 

■ 1800 UTC 17 January, 

■ 0000 UTC 18 January, 

• 0600 UTC 18 January, 

■ 1200 UTC 18 January, 

The 12-hr forecast grid (0000 UTC 18 January) used for boundary conditions was found to contain bad 
data that caused RAMS to stop running at the 7-h forecast time (1900 UTC), To fix the problem, we substituted 
the 24-h Eta forecast grid from the previous forecast cycle. RAMS ran to completion widi these data. 

Other meteorological data were available to characterize the weather conditions at the time of the accident. 
This data included the following: 

• WINDS tower data (available at 5-minute intervals) 

• Surface and buoy observations (available hourly) 

• Rawinsonde data (available at 1 128, 1243, 1453, 1548, 1613, and 1815 UTC) 

• 50 MHz profiler data (available hourly) 

• 915 MHz profiler data (including Radio Acoustic Sounding Systems virtual temperature data) 
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2.2.3. Satellite data 

Satellite data from the GOES-8-east satellite data was available for 17 January 1997. The satellite photos 
present on this day. The satellite data also used to show the location of the 

explosion clouds as they advected downwmd. 

2.3. Model description/configuration 

2J.1. ERDAS/PROWESS development 

ERDAS was developed by Mission Reseamh Consomdon (MRCRASTER dm»n 
pro, owe system ,0 provide 7 “ 8 en=yje^o^ 

„U :rr^.'Sc^:pr„ mode,™ 

“ L“ e *ob addressing before die ERDAS oonid become 

Operational (Evans 1996). 

Ti, Pcrallpliyed RAMS Operational Weather Simulation System (PROWESS) was developed by MRC 
during L (Tron*»k e, ai. 1996). The goafc of the PROV^SS proi». were m 

. ^ F K tk nf thp narallelized RAMS code to run on a cluster of low-cost workstations. The model 

rd.?“deX^.--cdo« a^ 

The nex. generation of E^AS w„ «HosUd imo 

Support Replacement (MARSS-REPL) system. MAR c nn nfMARSS and a prototype of the MARSS- 

primary toxic/hazard diffusion prediction sys^m. An earlier version of MARSSm^^^^ ^ 

diffiision moLl Ocean EreezeCIry Oolch. Users can abo display vartotis meutorologtcal displays. 

Some of die prnnary features of dte upgnuied «td certified system (called ERDAS-REPL), deliyeted in late 
1998, are; 

. An upgraded LAN using a higher speed distribution mechanism known as Fiber Distributed Data 
Interface to support the larger data transfers. 

• A series of Hewlett-Packard K460 workstations used as dedicated RAMS processors utilizmg a 
Symmetric Multi-Processor architecture. 

• RAMS with; 

■ Full microphysics activated, 

- Four nested grids with the largest grid covering the southeast United States, and 
. 1 ,25-km grid spacing on the fine grid covering approximately 1 00 km by 80 km centered 

on CCAFS/KSC. 

. Additional meteorological data sources now include 50 MHz Doppler rad^ wind Profiler, 9 15 
MHz radar wind profilers, rawinsondes, coastal buoys, surface stations, wind towers, and gndded 
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National Centers for Environmental Prediction (NCEP) ETA model grids. New data quality 
control procedures were implemented for observational data. 

• The system will now run Rocket Exhaust Effluent Diffusion Model (REEDM) and BLASTX 
models, along with Ocean Breeze/Dry Gulch (OB/DG), HYPACT/RAMS. 

23.2. RAMS-ERDAS 

The AMU conducted the Delta II modeling analysis by running RAMS for 17 January 1997 and using the 
forecast meteorological data to drive the diffusion model HYP ACT. The RAMS-ERDAS configuration was the 
same as the configuration used for the daily operation of RAMS in the prototype ERDAS. This configuration 
has been set since the prototype ERDAS was installed in the ROCC in 1994. The key points of this 
configuration are: 

• RAMS version 3a 

• Microphysics inactive (NLEVEL=1) 

• 3 nested grids 

■ Coarse grid: 60 km spacing, 2220 x 2 100 km domain 

■ Medium grid: 15 km spacing, 495 x 555 km domain 

■ Fine grid: 3 km spacing, 108 x 108 km domain 

• Vertical grid with 10-m lowest grid point on fine grid and expanding in depth upward 

• Twice-daily RAMS runs initialized at 0000 and 1200 UTC with hourly forecast output 

• Input data used to initialize the model are obtained from MIDDS and include: 

■ NCEP ETA data 

■ rawinsondes 

■ surface and buoy data 

• local CC AFS/KSC WINDS tower data 

A complete listing of the file RAMSIN containing the input parameters for running RAMS is provided in 
Appendix A. 


2.33. RAMS-PROWESS 


The RAMS-PROWESS runs were made using the version of RAMS on the PROWESS workstations. The 
PROWESS workstations consist of one IBM RS/6000-370 and seven IBM RS/6000-250 workstations which run 
the parallel version of RAMS. The key points of this configuration are: 

• RAMS version 4a 

• Microphysics active (NLEVEL=3) 

• 4 nested grids 

■ Coarse grid: 72 km spacing, 2376 x 2088 km domain 

■ Medium grid: 1 8 km spacing, 594 x 666 km domain 

■ Fine grid: 6 km spacing, 222 x 222 km domain 

■ Finest grid: 1.5 km spacing, 61.5 x 85.5 km domain 

• Vertical grid with 38-m lowest grid point 

• Twice-daily RAMS runs initialized at 0000 and 1200 UTC with hourly forecast output 
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• Hourly output beginning at initialization time of 1 200 UTC 

• Input data used to initialize the model is obtained from MIDDS and includes 

• NCEP ETA data 

■ Rawinsondes 

■ Surface and buoy data 

■ Local CCAFS/KSC WINDS tower data 

A complete listing of the file RAMSIN contaming the input parameters for running RAMS is provided in 
Appendix B. 


2.3.4. HYPACT-ERDAS and HYP ACT-PROWESS 

The primary model used for computing dispersion estiinates in ^R^AS is 
advanced Lagrangian particle dispersion model. Dispersion m the Lagrangian mode ^ 

^alSng a t^"et of partmles- Subsequent pos.tions of each particle are computed from the relatmn. 

X[t+At] = X[t] + [u + u'] At 
Y[t +At] = Y[t] + [v + v ] At 
Z[t +At] = Z[t] + [w + w' + Wp] At 

where u v and w are the resolvable scale wind components which are derived from Ae ^S ^ field 
and u'^v', a^d w' are the random subgrid turbulent wind components deuced from RAMS. The Wp tenn 
terminal velocity resulting from external forces such as gravitational settlmg. 

REEDM predicts plume rise and downwind ^ 

ERDAS, REEDM produces the source term which is used by HYP AC p P 

resulting downwind concentrations. 

For modelmg tstagTe Ss^^eirc^rJ^f fi^dT ^d 

HYPACT can model any number of sources which are specified anywhere in the domain and configured ^ 
pomX a^eror The emissions from these sources can be instantaneous, mtermritent, or 

continuous and the pollutants can be treated as gases or aerosols. 

To simulate the Delta II explosion plume using HYPACT required that we modify the HYPACT input 
JT ' A within the FRDAS dispersion function. For launch plume dispersion scenario , 

of « ^ vertical colurtrt,. repceaenia «.= »»,» of an en.i.n=d k,nnch plunae p^ducf REEDM 
computes the source term based on a number of factors such as. 

• type of launch (e.g. nominal, abort conflagration, or abort deflagration) 

. type of vehicle (e.g. Atlas, Delta II, Titan IV, Shuttle, etc.) 

• vertical stability (i.e. stabilization height). 


2.3.5. REEDM 

Prior to each land-based launch, Range Safety, in conjunction with their support contractor, ACTA Inc 

Zrtir, 

fcT^LTf SroM Sme" S^Sefy'ta cont.nued B upennie REEDM and uses d.e toes, 

version available to evaluate launch cloud dispersion. 
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Toxic substances are identified on the EPA list of Extremely Hazardous Substances (EHS). The Delta II 
contained or produced the following EHS chemicals as propellants or products of solid rocket propellant 
combustion. 

• Hydrazine 

• Unsymmetrical dimethyl hydrazine (UDMH) 

• Nitrogen tetroxide 

• Hydrogen chloride 

Aerozine-50 is the Delta II second stage hypergolic fuel and is a 50-50 mixture of hydrazine and UDMH. 
Nitrogen tetroxide is the Delta II second stage hypergolic oxidizer. The reaction of the solid propellant involves 
aluminum oxide as the oxidizer and ammonium perchlorate as the fuel to produce hydrogen chloride along with 
other non-EHS reaction products. 

REEDM models the toxic cloud formation process in four significant steps: 

• cloud formation 

• cloud rise (buoyant cloud rise) 

• cloud stabilization (at neutral buoyancy after heat exchange with the atmosphere) 

• cloud transport and dispersion 

A single rawinsonde sounding provides the input meteorological data for REEDM, which incorporates vertical 
but not horizontal wind shear. Thus the sounding is assumed to represent the flow across the entire area-wide model 
domain. 

REEDM predicts the ground trace or net ground level concentration of hazardous material arising due to its 
dispersion from the elevated cloud. This ground-level isopleth does not portray the cloud’s passage in terms of 
time, but lists range and bearing from the source, peak concentration at that range and bearing, and, cloud arrival 
and departure times at that range and bearing. REEDM results for the Delta II scenario are presented in Section 
3.2.1. 
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3. Results 

The results of the Delta H explosion plume analysis are presented in this section. The observati^ which 
provided the verification and which were used to compare with modeled data are presented fet. T^e 
Lluded meteorological, visual (video tape) observations, sateUite 

used in this analysis are presented next. These data included the predictions produced by REEDM, RAMS, and 
HYPACT. 

3.1. Observed data 


3.1.1. Meteorology on 17 January 1997 

A cold front had passed through the Cape Canaveral area approxi^tely 24 houK prior to the Delta II 
launch and had pushed through all of Florida, most of the Bahamas, and through Ac eastern 
(Figure 5). Therefore, Ae prevailmg weaAer conAtions at Cape Canaveral on H January were^st^onta 
na^. Across Ae Florida peninsula, Ae winds at Ae surface were generally ^m Ae norA and 
Ae cold dense air moved souAward. A high pressure area centered over Ae Louisima area w^build g 
eastward mto Florida. The colder air mass was relatively shaUow at launch tune as subsiding aur behmd Ae 
front formed a well-defined temperature Aversion near 900 m. Skies were clear over Ae launch site while 
scattered to broken stratocumulus clouds were located over Ae ocean further east. Low temperat^s m Ae area 
Aat mommg were A Ae upper 30’s. The Shuttle LandAg FacUity reported a low temperature of 36 F. 

A graph showAg Ae potential tempetature profile from Ae Arec Cape C^v^ rawAsondes on 17 
January are presented A Figure 6. The graphs show Ae very strong ™ 

meters at 1128 UTC. Also shown are Ae potential temperature profiles from 1613 and 1815 mc. At Aese 
times, Ae Aversion base had risen to approxAiately 900 meters as Ae surfitee warmed 5 to J tC A ^h 
showAg Ae wAd speed and wAd direction at 1815 UTC show Ae inark^ difference . » to i n m c ' 

below Ae Aversion (Figure 7). The northwesterly wAds below Ae mversion were fauly strong at 8 to 
but Acreased A speed and Aen backed to become more westerly above Ae Aversion. 

Table 1 presents Ae hourly wAd speed and direction data at 150-meter Acrements from 2000 to 5000 m as 
collect^ To MHz profiler at KSC. The profiler data AAcate Ae wAds at the 2000-meter^^^^^^^^ 

n<mhwesUTly with speeds of 13 » 18 knots « 1600 UTC but tb« beetme ^ wti*, *^S*e 

during Ae next five hours. The wAds above 2700 m were westerly wiA speeds of 20 to 30 knots durmg Ae 

period 1600 to 2100 UTC. 
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Table I . Wind speed and direction from the 50-MHz profiler at KSC on 17 January 1997. 


Time 

(UTC) 

1600 

1700 

1800 

1900 

2000 

2100 

Height 


Direct. 

(deg) 

Speed 

(^^) 

Direct. 

(deg) 

Speed 

(1^) 

Direct. 

(deg) 

Speed 

(kts) 

Direct. 

(deg) 

Speed 

(kts) 

Direct. 

(deg) 

Speed 

(lets) 

Direct. 

(deg) 

Speed 

(kts) 

2011.0 

318 

13.5 

268 

20.2 

272 

20 

284 

18.7 

285 

17.4 

290 

16.8 

2161.0 

316 

15.3 

276 

20.6 

279 

19.5 

286 

19.2 

269 

18.3 

288 

15.6 

2311.0 

288 

18.6 

289 

16.8 

281 

21.0 

284 

18.7 

284 

17.9 

285 

16.5 

2461.0 

287 

20.4 

288 

18.2 

288 

17.8 

284 

17.4 

282 

17.6 

290 

17.9 

2611.0 

286 

19.7 

285 

18.6 

285 

18.2 

282 

18.2 

281 

18.4 

292 

19.0 

2761.0 

279 

20.0 

278 

20.5 

283 

20.9 

285 

19.9 

285 

19.8 

291 

19.6 

2911.0 

274 

20.8 

273 

20.4 

282 

21.3 

283 

20.5 

283 

21.4 

284 

21.9 

3061.0 

278 

20.5 

270 

20.5 

271 

21.1 

275 

21.2 

276 

22.7 

289 

21.9 

3211.0 

277 

20.9 

270 

20.7 

270 

21.7 

274 

21.5 

277 

21.6 

284 

23.0 

3361.0 

276 

21.7 

275 

21.2 

277 

22.7 

278 

21.8 

277 

22.1 

281 

23.1 

3511.0 

277 

23.1 

278 

22.4 

281 

23.8 

278 

23.5 

278 

23.1 

277 

23.7 

3661.0 

280 

24.3 

283 

25.1 

282 

26.7 

278 

25.2 

276 

24.8 

276 

25.1 

3811.0 

281 

25.4 

284 

26.7 

283 

27.1 

279 

26.0 

275 

25.7 

275 

25.7 

3961.0 

276 

26.0 

280 

26.7 

282 

26.7 

276 

26.1 

272 

27.1 

271 

26.5 

4111.0 

268 

25.9 

270 

26.2 

272 

27.3 

268 

26.7 

267 

27.6 

268 

28.3 

4261.0 

266 

26.0 

264 

26.5 

266 

28.5 

266 

27.8 

264 

28.2 

262 

28.1 

4411.0 

267 

25.9 

263 

27.1 

263 

27.6 

262 

27.1 

263 

28.6 

260 

28.4 

4561.0 

263 

24.9 

259 

25.3 

260 

26.8 

259 

26.6 

260 

28.1 

259 

29.0 

4711.0 

261 

24.7 

258 

25.7 

259 

27.0 

258 

27.2 

256 

28.9 

262 

29.7 

4861.0 

260 

25.6 

257 

26.4 

258 

28.5 

256 

27.9 

260 

30.5 

264 

30.6 

5011.0 

262 

27.4 

263 

28.7 

261 

29.7 

261 

29.8 

264 

30.2 

263 

32.7 


3.1.2. Visual observations 

The Delta II explosion was captured on videotape by numerous Air Force and amateur video cameras. 
Most of the Air Force video focused very closely on the exploding rocket since these pictures were used for the 
investigation to determine the cause of the explosion. The amateur video provided a better look at the clouds 
resulting from the explosion since they were taken from a distance and showed a full panorama of the upper 
and lower clouds. Pictures extracted from on amateur video are presented in Figures 8-10. The person who 
shot the video and the exact location are unknown. However, the video was taken from the beach in Cocoa 
Beach. Bionetics was responsible for archiving all amateur and Air Force video of the Delta II explosion. 

The video shows the huge lower cloud that was produced when the nine solid propellant motors and first 
stage exploded at approximately 13 seconds. The video also shows the upward movement and subsequent 
destruction/explosion of the Delta II’s stage 2 and 3. The pictures show the upper plume to be darker and much 
smaller than the lower plume. The captured frames in Figures 8 to 10 show the two distinct clouds very well 
along with some other notable features. The inversion that existed at 900 meters was clearly evident as it 
capped the lower cloud and allowed very little of the buoyant smoke and gas to rise past the inversion. 

One very key feature shown in the video and in the captured figures was the small dark orange cloud at the 
top of the lower cloud (Figure 10). This dark orange cloud indicates the presence of nitrogen tetroxide (N2O4) 
that was emitted when the second stage was destructed. Ha and Deane (1998) indicated that Stage 2 contained 
8759 Ibm of N2O4. Since this cloud was at the top of the lower cloud, the video indicates that this cloud was not 
completely trapped below the inversion, as was most of the lower cloud. This observation is important in the 
determination of the transport of the N2O4 because of its toxic nature (Sittig 1985). The following sections 
describe in detail the transport of the lower cloud as tracked by satellite and radar but the fate of the small 
orange cloud at the top of the lower cloud cannot be directly tied to the large lower cloud. The bulk of the N2O4 
was likely contained in the upper cloud and was transported toward the east-southeast. 
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Figure 8. The two clouds produced from the explosion of the vehicle a few seconds after the explosion on 17 
January 1997. This picture was captured from amateur video taken along Cocoa Beach. 
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Figure 9. The two explosion clouds at approximately 16 seconds after the explosion. Note the orange area 
located at the top of Plume 1 . 



Figure 10. Zoomed in view of explosion clouds at approximately 12 seconds after the explosion. The small 
orange cloud at the top of the lower cloud most likely contains nitrogen tetroxide. 
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3.U. Satellite data 


GOES satellite data for two times following die explosion is shown in Figures 1 1 md 12 (Brandli 1997). 
These satellite pictures show the lower explosion cloud as it moved south over the Melboime ar^. Figi^e 11 
shows the cloud as it reaches the coastline in Cocoa Beach at 1645 UTC while Figure 12 shows the cloud as it 
reaches the Melbourne area at 1715 UTC. 



Figure 11. GOES satellite image of Delta II explosion cloud at 1645 UTC 17 January 1997. (from Brandli 
1997). 
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Figure 12. GOES satellite image of Delta II explosion cloud at 1715 UTC 17 January 1997. (from Brandli 
1997). 


3.1.4. Radar data 

For the 1 7 January case, the WSR-88D clearly displayed the echoes resulting frx)m the rocket explosion. As 
discussed earlier, the display showed two distinctly separate plumes. Information on the two plumes was 
obtained for each 10-minute scan by using the 88Display information query function. This information is 
presented in Tables 2 and 3. Much of the data contained in the table was obtained by manually analyzing each 
time and elevation scan using the information query function. The scan start time, stop time, and elevation 
angle were displayed for each scan. We zoomed in on the area of interest to show the explosion clouds and 
found the location of the maximum reflectivity using the color contouring. By selecting this point with the 
mouse, we obtained the vital information such as location and height for that time and elevation. The actual 
maximum reflectivity value was estimated fi“om the color contours. The plume dimensions such as downwind 
distance and alongwind and crosswind length were estimated using the cross-section tool of 88Display. 

The plume top and bottom were found using the information queiy function for each elevation scan and 
finding the highest and lowest point of the cloud of interest. It should be noted that the heights of the clouds 
determined from this technique have a large uncertainty for reasons discussed earlier. 

The Delta II abort cloud was first detected by the WSR-88D during the scan at 1635-1646 UTC, Figures 
13 and 14 show the composite reflectivity and cross section of the radar image at that time. The cross section 
shows how the plume split into two distinct clouds with an upper cloud to the east and a lower cloud to the 
south of LC-17. Figures 15 and 16 show the clouds two hours later at the scan ending at 1834 UTC. The upper 
and lower clouds were still clearly detected by the radar. The lower cloud was more distinguishable, located 
south of the Melbourne area. 

The composite reflectivity figures for each of the scans arc shown in Figures 17 to 42. 


18 



Table 2. Data obtained from WSR-88D Level II archive for 17 January 1997 using the 88Display software. 
Data in this table show Plume 1 (lower cloud) information. 


Scan Scan Scan Lat. of Lon.of Height Max Plume 1 Plume I Plume 1 Plume 1 Plume I Plume I 
Time Time Time Peak Peak of Peak Reflect Center Center Down- Direction Along Cross- 

Start Start Stop Ut Lon. wind fromLC- wind wind 

Dist 17 length Dist. 

(hhmmss) (hrs) (hhmmss (dcgN) (degW) (km) (dBZ) (degN) (degW) (km) (deg) (km) (km) 

) 

162552 16.4456 163639 28.4355 80.5641 0.444 45 28.4319 80.5641 1 180 2 3 

163541 16.5947 164628 28.3839 80.5692 0.37 53 28.3863 80.5678 7 185 3 5.5 

164531 16.7586 165616 28.3212 80.5733 0.29 47 28.3212 80.5802 13 190 5 1 1 

165519 16.9219 170604 28.2667 80.5763 0.22 45 28.2477 80.6076 21 190 8 9 

170538 17.0939 171519 28.2109 80.5748 0.17 39 28.1708 80.6300 30 190 11 13 

171527 17.2575 172507 28.1554 80.5708 0.16 33 28.1589 80.5788 33 190 10 11 

172515 17.4208 173455 28.0677 80.6284 0.15 33 28.0849 80.5927 42 200 10 25 

173503 17.5842 174443 28.0356 80.6126 0.29 31 28.0303 80.5707 48 200 12 28 

174451 17.7475 175431 27.9897 80.5288 0.18 28 27.9844 80.5528 51 180 11 35 

175440 17.9111 180419 27.9726 80.5509 0.21 34 27.9612 80.5249 59 180 11 35 

180428 18.0744 181407 27.8970 80.4882 0.33 38 27.8970 80.4882 62 175 15 38 

181416 18.2378 182356 27.8581 80.5012 0.34 30 27.8581 80.5012 66 170 15 37 

182404 18.4011 183344 27.7985 80.4780 0.46 28 27.7985 80.4780 72 170 15 47 

183352 18.5644 184332 27.8166 80.5704 0.39 28 27.7707 80.4855 75 180 14 45 

184341 18.7281 185320 27.7805 80.5778 0.44 27 27.7380 80.5003 77 180 14 52 

185329 18.8914 190309 27.7454 80.5797 0.49 25 27.6865 80.4763 84 180 12 40 

190317 19.0547 191258 27.7159 80.5945 0.532 23 27.6472 80.4838 90 180 15 55 

191306 19.2183 192246 27.6865 80.5871 0.579 21 27.6177 80.4802 95 180 16 60 

192254 19.3817 193234 27.6505 80.6056 1.45 18 27.5981 80.5134 95 180 14 65 

193242 19.5450 194222 27.6014 80.5872 1.71 18 27.5948 80.5282 97 180 15 47 

194231 19.7086 195316 27.5785 80.5798 0.76 17 27.5785 80.5798 97 180 16 45 

195219 19.8719 200305 27.5457 80.5872 0.81 17 27.5457 80.5872 101 180 20 45 

200207 20.0353 201253 27.5359 80.6242 0.83 17 27.5359 80.6242 102 185 18 50 

201156 20.1989 202241 27.5065 80.6057 0.88 16 27.5065 80.6057 105 190 22 52 

202144 20.3622 203230 27.4901 80.6352 0.91 14 27.4901 80.6352 110 190 20 45 


Plume 1 Plume I 
Top Bottom 


(km) (km) 

~ 097 

2.1 0.35 

1.9 0.27 

1.2 0.2 

1.1 0.08 

1.1 0.04 

1.1 0.04 

0.91 0.06 

0.99 0.09 

1.13 0.13 

1 .09 0. 1 7 

1 0.2 

1.08 0.24 

1.22 0.3 

1.35 0.31 

1.57 0.37 

1.52 0.41 

1.51 0.45 

1.35 0.45 

1.13 0.56 

1.15 0.54 

1.18 0.59 

1.34 0.59 

1.38 0.69 

1.31 0.69 


Table 3. Data obtained from WSR-88D Level II archive for 17 January 1997 using the 88Display software. 
Data in this table show Plume 2 (upper cloud) information. 


Scan Time 
Start 

(hhmmss) 

Scan 

Time 

Start 

(hrs) 

Scan 

Time Stop 
(hhmmss) 

162552 

16.4456 

163639 

163541 

16.5947 

164628 

164531 

16.7586 

165616 

165519 

16.9219 

170604 

170538 

17.0939 

171519 

171527 

17.2575 

172507 

172515 

17.4208 

173455 

173503 

17.5842 

174443 

174451 

17.7475 

175431 

175440 

17.9111 

180419 

180428 

18.0744 

181407 

181416 

18.2378 

182356 

182404 

18.4011 

183344 

183352 

18.5644 

184332 

184341 

18.7281 

185320 

185329 

18.8914 

190309 

190317 

19.0547 

191258 

191306 

19.2183 

192246 

192254 

19.3817 

193234 

193242 

19.5450 

194222 

194231 

19.7086 

195316 

195219 

19.8719 

200305 

200207 

20.0353 

201253 

201156 

20.1989 

202241 


Plume 2 
Center Lat. 


Plume 2 
Center Lon. 


(degW) 

80.5011 

80.4089 

80.3391 

80.2657 

80.1466 

80.0690 

80.0824 

80.0044 

79.9306 

79,8300 

79.7757 

79,6951 

79.6354 

79,4792 

79.4477 

79.3665 

79,2517 

79.1706 

79.1867 

79.2281 

79.7336 

79.6710 

79.6310 

79.6315 


Plume 2 
Down-wind 
Distance 


Plume 2 
Direction - 
from 
LC-17 


110 

no 

no 

115 

115 

120 

120 

115 

no 

no 

no 

no 

115 

120 

120 

120 

120 

120 

120 

125 

135 

135 

135 

135 


Plume 2 
Cross-wind 
Distance 
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Figure 13. Composite radar image from Melbourne WSR-88D for 10-minute scan beginning at 1635 UTC. A-B 
line indicates cross section shown in Figure 14. 



Figure 14. Cross section of composite radar image for 1635 UTC. East-west cross section location shown in 
Figure 13. 
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Figure 15. Composite radar image from Melbourne WSR-88D for 10-minute scan beginning at 1833 UTC. A-B 
line indicates cross section shown in Figiue 16. 



Figure 16. Cross section of composite radar image for 1833 UTC. East-west cross section location shown in 
Figure 15. 
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3.I.4.I. 


Cloud 1 


The lower cloud, designated Cloud 1, moved to the south over the next four hours after the explosion as 
the WSR-88D radar clearly tracked the smoke and particles contained in the cloud. The reflectivity reached a 
maximum of 53 dBZ in the 1635-1646 UTC scan (Figures 13 and 19) but then gradually decreased to 
approximately 14 dBZ four hours later in the 2021-2032 UTC scan (Figure 42). The shape of the cloud as 
shown by radar appeared rounded initially but then became oval-shaped in the crosswind direction at 
approximately one hour after the explosion. The cloud maintained this oval shape for the next three hours as it 
became elongated in a northwest-southeast orientation. The analysis of the 1635 UTC scan (Figures 13 and 19) 
indicated the low-level cloud was approximately 3-km long, 2-km wide, at least 2-km thick and was located 1 
km from LC-17. The cloud thickness is likely greatly exaggerated because of radar beam width and 
superrefraction effects. The direction the cloud moved was initially south and it maintained this southerly 
movement throughout the four hours. The cloud moved to intersect the coastline at approximately 1650 UTC 
(Figure 20) as it moved over the beaches between Cocoa Beach and Patrick AFB. It continued to move south 
over the Melbourne area and was centered over the Melbourne airport at approximately 1735 UTC (Figure 25). 
The cloud moved from Brevard County into Indian River County at {^)proximately 1830 UTC (Figures 15 and 
30). 

Ground fires produced smoke plumes which were detected by flie WSR-88D on 17 January. These fires as 
observed on the radar displays give a good indication of the surfau^ level winds. The display at 1804 UTC 
(Figure 28) shows the smoke plume produced by the fires over Cs^ Canaveral that were started from the 
falling solid rocket fuel. The low level winds affecting this plume were fix>m the north-northwest. It is 
important to note that this low-level smoke plume stays offshore and does not move westward as did the Delta 
II explosion cloud which was aloft. The plume from another fire (unrelated to the Delta II explosion) located in 
western Indian River County (southwest of the lower Delta II cloud) moves south and south-southwest. The 
low-level winds in this area were from the north and north-northeast. 

The height of the cloud top and bottom as determined by the 88Display program on the Level II data is 
shown in Table 2. As discussed earlier these values are not very accurate because of errors caused by refiuction 
of the radar beam through a strong temperature inversion and because the cloud passed directly over the 
Melbourne radar site. However, the cloud height data does provide some useful information. The plume top 
was initially determined to be 2700 meters which is obviously too high as determined by the capping inversion 
characteristics. The estimated plume top dropped to 1200 meters at 1655 UTC which more closely matches the 
other observations. The reason for the large initial discrepancy is probably due to the beam width and the beam 
inversion problem. The cloud bottom was initially measured at 970 meters and dropped to a low of 40 meters at 
1725 UTC as the cloud passed over the radar site and then rose to 690 meters at 201 1 UTC. The change over 
time of the height of the plume top and bottom most likely reflects the radar beam angle since the heights are 
directly proportional to the distance from the radar. A graph showing the plume height measurements is 
presented in Figure 43. 
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3.1A2. Cloud 2 


The Delta II stages 2 and 3 along with the payload survived the initial explosion and drifted upward to just 
over 750 meters where destruct signals were sent at 22.4 seconds. This formed a second cloud which broke 
through the inversion. The upper cloud moved to the east and southeast over the four hours after the explosion 
at an average speed of approximately 13 m s ^ This cloud was detached from the lower cloud and moved at a 
different speed and direction. The radar detected peak reflectivities of approximately 34 dBZ at 1645 UTC 
(Figure 20). The reflectivities slowly decreased over time, then dropped below 30 dBZ at 1715 UTC (Figure 23) 
and then below 20 dBZ at 1754 UTC (Figure 27). The cloud was moved by the strong westerly winds in the 
layer above the inversion between 1200 and 2000 meters. The radar showed that the cloud had a somewhat 
undulating shape as it was stretched from west to east. 

The height of the cloud top and bottom as determined by the 88Display program on the Level II data is 
shown in Table 3. The initial scan of the upper level cloud at 1625 UTC indicated the cloud was approximately 
4-km long, 4-km wide, 1.5-km thick and was located 6 km east of LC-17. The errors associated with the 
refraction of the radar beam through the strong inversion that caused estimation errors for the height of the 
lower cloud were more prominent in determining the height of the upper cloud. The radar indicated that the 
initial vertical plume dimensions ranged from a cloud bottom of 1700 meters up to a cloud top of 3200 meters 
(Figure 14). During the next three hours the radar detected that the cloud top rose to a height of 4400 meters at 
1833 UTC while the cloud bottom lowered to a height of 450 meters at 1804 UTC (Table 3). The 
measurements indicate a maximum depth of the upper cloud of over 3700 meters at 1833 UTC. Visual 
observations discussed in Section 3.1.2 did not support these radar depth determinations. Therefore, we 
seriously question this cloud’s radar height measurements which are shown in a graph in Figure 44. 
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Figure 17. WSR-88D composite scan beginning at 
1616 UTC. 


Figure 18. WSR-88D composite scan beginning at 
1625 UTC. 



Figure 19. WSR-88D composite scan beginning at 
1635 UTC. 


Figure 20. WSR-88D composite scan beginning at 
1645 UTC. 






Figure 21. WSR-88D composite scan beginning at 
1655 UTC. 


Figure 22. WSR-88D composite scan begiiming at 
1705 UTC. 




Figure 23. WSR-88D composite scan beginning at 
1715 UTC. 


Figure 24. WSR-88D composite scan beginning at 
1725 UTC. 






Figure 25. WSR-88D composite scan beginning at 
1735 UTC. 



Figure 27. WSR-88D composite scan begirming at 
1754 UTC. 


Figure 26. WSR-88D composite scan beginning at 
1745 UTC. 



Figure 28. WSR-88D composite scan beginning at 
1804 UTC. 






Figure 29. WSR-88D composite scan beginning at 
1814 UTC. 


Figure 30. WSR-88D composite scan begiiming at 
1824 UTC. 



Figure 3 1 . WSR-88D composite scan beginning at 
1833 UTC. 


Figure 32. WSR-88D composite scan beginning at 
1843 UTC. 







Figure 33. WSR-88D composite scan beginning at 
1853 UTC. 



Figure 35. WSR-88D composite scan beginning at 
1913 UTC. 



Figure 34. WSR-88D composite scan beginning at 
1903 UTC. 



Figure 36. WSR-88D composite scan begiiming at 
1922 UTC. 
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Figure 37. WSR-88D composite scan beginning at Figure 38. WSR-88D composite scan beginning at 

1932 UTC. 1942 UTC. 



Figure 39. WSR-88D composite scan beginning at Figure 40. WSR-88D composite scan beginning at 

1952 UTC. 2002 UTC. 
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Figure 41. WSR-88D composite scan beginning at 
2011 UTC. 



Figure 42. WSR-88D composite scan beginning at 
2021 UTC. 
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Figure 43. The plume height of the lower cloud (Plume 1) over the 4-h period. Heights were determined using 
the 88Display software tool. Note that the inaccuracies of these height determinations is high due to radar beam 
width spread and radar beam path through the strong inversion. However, the determination of error bars was 
not possible using the 88DispIay tool. 



Figure 44. The plume height of the upper cloud (Plume 2) over the 4-h period. Heights were determined using 
the 88Display software tool. Note that the inaccuracies of these height determinations is high due to radar beam 
width spread and radar beam path through the strong inversion. However, the determination of error bars was 
not possible using the 88Display tool. 
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3.2. Model data 


One of the primary goals of this study was to compare the model results with observations for 17 January 
1997, This section presents the model data from RAMS, HYP ACT, and REEDM for this day. Comparison with 
observational data discussed previously is provided. 

3.2.1. REEDM results 

These REEDM results were obtained from ACTA (Figures 45 and 46). A rawinsonde released at 1613 UTC 
provided data for input to REEDM. The REEDM vertical profiles of temperature, wind speed, wind direction, 
and potential temperature are shown in Figure 46, Winds were northerly with wind speeds near 1 1 m s'^ below 
the inversion, becoming northwesterly then westerly and stronger above the inversion. 



Figure 45. REEDM plots from the 1613 UTC rawinsonde just prior to the Delta II rocket failure on 17 January 
1997. The REEDM run was for a conflagration at T+15 seconds simulating transport and ground-level 
concentration of HCl gas from burning solid rocket motors. 


Following the Delta El accident, REEDM modeled only the lower toxic cloud, which was a potential health 
threat at the surface. The sounding detected northerly or northwesterly winds below the inversion, so the REEDM 
calculation drove the cloud southward, with landfall south of Melbourne. The peak HCl concentration predicted by 
REEDM as the plume moved southward was on the order of 1.0 ppm for the first 30 min; the highest peak (1.27 
ppm) was predicted to occur 1.5 km downwind from the pad. 

At the right of the Figure 46 is the predicted profile of the solid rocket propellant cloud as viewed from a 
position west of LC-17. Winds above and below the inversion sheared the cloud, moving the portion below the 
inversion toward the south, and the portion above the inversion toward the southeast. The REEDM predicted 
plume track and 1-ppm contour of ground-level HCl concentrations are shown in Figure 45. Not surprisingly, 
given the meteorological conditions shown in Figure 46, the plume was predicted to move almost due south 
over the ocean and not make landfall until south of the Melbourne area. 
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The results presented here were from REEDM 7.07 which was the version in use at the time of the Delta II 
explosion. Newer versions of REEDM (e.g. REEDM 7.08 and REEDM 7.09) produce different results and 
different ground level concentrations. Range Safety has continued to upgrade REEDM as the result of new 
scientific findings and uses the latest version available to evaluate launch cloud dispersion (Nyman 1999). 



Figure 46. REEDM cross-section plot showing rawinsonde data from 1613 UTC on the left and REEDM 
source term predictions on the right. 

3.2.2. RAMS results from ERDAS and PROWESS 

The description and configuration of the two versions of RAMS used in this study were provided in 
Section 2.3. Discussion and figures showing the results of those RAMS runs are presented here. 

The RAMS runs were initialized at 1200 UTC, approximately 4.5 hours before the Delta II explosion. The 
meteorology during this day did not change significantly because of the presence of the post-frontal regime 
with the continued cold air advection. There were no sea breezes, river breezes, thunderstorms, or other 
significant local-scale phenomena for RAMS to forecast for the Cape Canaveral area. 

Graphs comparing ERDAS RAMS wind data with observed data are presented in Figures 47 and 48. 
These figures show the observed vertical wind profiles at 1613 UTC and the predicted vertical wind profiles at 
1600 UTC. RAMS predicted the wind speed profile accurately through the mixing layer to the inversion at 900 
meters and then accurately predicted the gradual increase to 25 m s ' at 4500 meters. RAMS predicted the wind 
direction more westerly than northwesterly in the layer from the surface to approximately 500 meters but then 
accurately predicted the wind shear at the 900-meter inversion and the westerly winds aloft. 

The potential temperature profile predicted by RAMS at 1600 UTC closely matched the observed profile at 
1613 UTC (Figure 49). The base of the predicted inversion was only slightly lower than the base of the 
observed inversion. 

Maps showing RAMS forecast data from ERDAS and PROWESS for the period 1500 UTC to 2000 UTC 
are shown in Figures 50 - 56. Data are presented that represent wind flow in Ae different layers important in 
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the plume analysis. The wind flow as shown by streamline analysis for different times and different levels are 

presented in Figures 50 to 53. The streamlines indicate the wind direction at a point in time and snace 

Wind Sf»ed(m-sec-1) ^ 

0 5 10 1$ 20 25 30 



Figure 47. Observed wind speed and direction from rawinsonde at 1613 UTC. Note that the observed wind 
directions from 0 - 15° are shown on the graph as greater than 360°. 

Wind Speed (m-sec-1) 

0 5 10 15 20 25 30 



Figure 48. Predicted wind speed and direction from ERDAS RAMS at 1600 UTC. Note that the observed wind 
directions from 0 - 1 5° are shown on the graph as greater than 360°. 
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Figure 49. Rawinsonde-observed potential temperature at 1613 UTC compared to model-predicted potential 
temperature at 1600 UTC. The model predictions were from ERDAS RAMS. 

The streamlines for the lowest RAMS level in both configurations indicate persistent northwesterly flow 
for the duration of the Delta II plume analysis over the Cape Canaveral region (Figs. 50, 51a and b). The 
models predicted a slight shift in flow from northwesterly to north-northwesterly from 1500 to 1900 UTC but 
the shift was not significant with respect to the plume transport in HYP ACT. 

The streamlines for the layers centered at 782 meters in ERDAS and 724 meters in PROWESS are shown 
in Figures 51c, 51d, and 52 for the Delta II analysis period. The 700-800 meter layer is important because this 
is the level just below the strong inversion that existed on 17 January. The winds in this level strongly 
influenced the transport of the lower plume. The analyses show that RAMS predicted northerly flow over Cape 
Canaveral at 1500 UTC. By 1700 UTC, RAMS predicted northeasterly wind flow in the region south of Cape 
Canaveral, over the ocean. The northeast flow became more pronounced at 1900 UTC. The ERDAS and 
PROWESS configurations generally agreed in predicting northeast winds south of Cape Canaveral but for the 
area north and west of Cape Canaveral, PROWESS RAMS predicted northerly flow while ERDAS RAMS 
predicted northeasterly flow across the entire Cape Canaveral region. 

For the layers well above the inversion at 1580 meters for ERDAS RAMS and 1699 meters for PROWESS 
RAMS, the streamline analysis from both models showed persistent west-northwesterly flow through the Delta 
II analysis period (Fig. 53). Both models showed similar wind flow at 1500 and 2000 UTC. The winds at this 
level influenced the transport of the upper cloud. 

The RAMS temperature predictions for the lowest levels in the ERDAS and PROWESS configurations are 
shown in Figure 54. The analyses are shown to compare the configuration differences and also to show the 
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change over time of the predicted low-level temperature structure. At 1500 UTC, both configurations show an 
east-west temperature gradient with colder temperatures over the land and warmer temperatures to the east over 
the ocean. It is difficult to compare the actual temperatures since the lowest level in ERDAS is centered at 10 
meters and the lowest level in PROWESS is centered at 35 meters. By 2000 UTC, the temperature gradient had 
shifted to north-south in both models and ERDAS RAMS, with its lowest level at 10 m, shows more influence 
of the land-water interfaces associated with rivers and islands. 


The vertical potential temperature structure of the atmosphere as predicted by RAMS is shown in Figure 
55. These figures show an east-west cross section of the potential temperature as predicted by the two 
configurations of RAMS at different times. The figures show that ERDAS RAMS predicted an elevated 
inversion at approximately 750 meters at 1500 UTC while PROWESS RAMS predicted the inversion at 
approximately 500 meters at 1500 UTC. By 2000 UTC the models lifted the inversion to approximately 1000 
meters for ERDAS RAMS and 750 meters for PROWESS RAMS. Rawinsonde data fi-om 1613 UTC indicated 
the base of the inversion was at approximately 900 meters (Figure 49). 
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Figure 50. S.reauiline forecast comparing ftom ERDAS and PROWESS on 17 Jannaiy 1997. Each 

figure is marked with its model configuration, time, and height. 
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a. ERDAS RAMS, 1900 UTC at 10 m. c. ERDAS RAMS, 1500 UTC at 782 m. 



b. PROWESS RAMS, 1900 UTC at 35 m. d. PROWESS RAMS, 1500 UTC at 724 m. 

Figure 51. Streamline forecasts comparing output from ERDAS and PROWESS on 17 January 1997. Each 
figure is marked with its model configuration, time, and height. 
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Figure 52. Streamline forecasts comparing output from ERDAS and PROWESS on 17 January 1997. Each 
figure is marked with its model configuration, time, and height 
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b. PROWESS RAMS, 1500 UTC at 1699 m. d. PROWESS RAMS, 2000 UTC at 1699 m. 

Figure 53. Streamline forecasts comparing output from ERDAS and PROWESS on 17 January 1997. Each 
figure is marked with its model configuration, time, and height. 
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ThJ^dark bold line in panels c and d shows the approximate location of the potential temperature cross sec 
shown in Figure 55 . 



35C0. 

3U00 

asoo 
2000 . 
« liCO 
1000 
500 
0 



■' " ■■■ — 

i 

) 

: 

i 

1 

I 

lS 



160 300. 220, 2^0 200 

X (klii) 

potent itil temp (K) 


y ^ -173. kiu 01/17 97 1500 UTC 


a. ERDAS RAMS, 1500 UTC. 
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Figure 55. Potential temperature (K) cross section forecasts comparing output from ERDAS and PROWESS on 
17 January 1997. Each figure is marked with its model configuration and time. The cross-sections intersect the 
Cape Canaveral Air Force Station along an east-west line as shown in Figure 54. The ERDAS cross-section is 
approximately 1 10 km wide and the PROWESS cross section is approximately 25 km wide. 
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3.2.3. HYP ACT results from ERDAS and PROWESS 


3.2.3.I. HYP ACT Configuration 

As discussed in Section 2.3.4, we modified the HYPACT input data configuration file to simulate the Delta 
II exoloSn^r^en we first ran REEDM within ERDAS for the Delta II explosion the somce term 
charTcterized as 21 different sources extending from 300 to 2100 m above the ground. Table 4 presents the 
source term data produced by REEDM for the 2 1 sources. 

Table 4. Location and dimensions of plume determined by REEDM and modeled by HYPACT. These plume 


Source 

latitude 

longitude 

height 

(m) 

x-size (m) 

y-size (m) 

z-size (m) 

shape 

1 

28.42615 

-80.55101 

349.9 

470.4 

470.4 

100.0 

ellipse 

2 

28.42622 

-80.55413 

449.9 

705.8 

705.8 

100.0 

ellipse 

3 

28.42408 

-80.55320 

508.3 

789.0 

789.0 

16.7 

ellipse 

4 

28.42531 

-80.55086 

524.9 

808.2 

808.2 

16.7 

ellipse 

5 

28.42695 

-80.54810 

550.0 

833.8 

833.8 

33.4 

ellipse 

6 

28.42695 

-80.54810 

583.5 

862.0 

862.0 

33.5 

ellipse 

7 

28.42695 

-80.54810 

614.5 

883.0 

883.0 

28.5 

ellipse 

8 

28.42695 

-80.54810 

643.0 

898.0 

898.0 

28.6 

ellipse 

9 

28.42695 

-80.54810 

671.5 

909.2 

909.2 

28.5 

ellipse 

10 

28.42695 

-80.54810 

742.9 

921.2 

921.2 

114.3 

ellipse 

1 1 

28.42695 

-80.54810 

856.3 

894.6 

894.6 

112.5 

ellipse 

12 

28.42695 

-80.54810 1 

925.1 

848.8 

848.8 

25.0 

ellipse 

13 

28.42695 

-80.54810 

950.1 ! 

826.0 

826.0 1 

25.0 ' 

ellipse 

1 A 

28 42695 

-80.54810 

975.1 

799.4 

799.4 

25.0 

ellipse 

15 

28.42695 

-80.54810 

1000.9 

767.4 

767.4 

26.7 

ellipse 

16 

28.42695 

-80.54810 

1028.6 

727.6 

727.6 

28.6 

ellipse 

17 

28.42695 

-80.54810 

1057.2 

679.4 

679.4 

28.6 

ellipse 

18 

28.42695 

-80.54810 

1085.8 

622.2 

622.2 

28.5 

ellipse 

19 

28.42695 

-80.54810 

1200.0 

399.8 

399.8 

200.0 

ellipse 

20 

28.42695 

-80.54810 

1450.1 

399.8 

399.8 

300.2 

ellipse 

21 

28.42695 

-80.54810 

2050.1 

399.8 

399.8 

899.8 

ellipse 


The visual and radar obseiwations of the plume mdicatea inai me piuu.c ..au uuv.v,-.., ...... •• - 

clouds To initialize more accurately HYPACT, we modified the source size and shape to more closely repose 

“sJr^eTs^rand shape of J actual plume. Therefore the ph^e w^ HY?ACT 

plume’s location and dimensions are presented in Table 5. A complete listing of the HYPACT configuration 

used to model the Delta II explosion is presented m Appendices C and D. 

Table 5. Location and dimensions of plume detennined ^om vis»l mdar ob^at™^ and modeled by 
... I HVPArT siTTiulatlon. 



BSSS 

longitude 

■BSH 

191 


^^9 



28.448 

-80.566 

600. 




ellipse 

2 

28.448 

-80.566 





ellipse 
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3.2.3J. 


HYPACT runs 


The HYPACT model was run for two scenarios. There were two reasons for the modeling exercise of 
comparing the observed and predicted plumes. The principal of the two reasons was to determine how well the 
modeled plume trajectories compared with the observed plume trajectories. The secondary reason for the 
exercise was to determine how the REEDM-predicted source term compared with the actual source term. To 
reduce the number of runs and the number of figures we displayed, we combined the trajectory and source term 
analysis in the ERDAS-HYPACT runs and we did not perform a source term analysis for the PROWESS- 
HYP ACT runs. Combining the runs did not hinder the results since we adequately assessed the trajectories of 
both runs and we were able to assess how the REEDM-predicted source term compeired with the actual source 
term. 

The figures presented in this report compare the predicted plume locations from ERDAS-HYPACT and 
PROWESS-HYP ACT with the observed plume locations from the WSR-88D radar. The comparisons presented 
are for each 10-minute HYPACT model time during the four hours that the plume was tracked. The HYPACT 
source terms were generated from the REEDM model using the REEDM function of ERDAS. REEDM 
generates the HYPACT source term for a rocket launch by creating a column that contains mass of the chemical 
species of interest. The column is composed of separate sources (volumes) of mass, which are 75 meters in 
thickness. We assessed the REEDM-predicted source term by observing how the modeled initial plume changed 
during the simulation and comparing its sh^e with observed plume. The ERDAS-HYPACT run was made 
leaving the mass in the entire column from the surface up to 2500 meters. We assessed the ERDAS-HYPACT 
trajectories by tracking only the lower and upper part of the plume and comparing those with the trajectory of 
the observed plume. 

The source term for the PROWESS configuration was modified slightly to simulate the observed source 
term since in our analysis we knew the initial shape of the plume. The plume was modified by removing the 
initial mass in the layers between 925 and 1550 meters. We assessed the PROWESS-HYP ACT trajectories by 
tracking the two plumes and comparing those with the observed plume. 

For the explosion time, both ERDAS-HYPACT and PROWESS-HYP ACT transported the upper and lower 
explosion clouds in different directions. The split transport occurred because the upper cloud was at a height 
where the winds were strong and from the west while the lower cloud was trapped below a strong inversion with 
winds were from the north. 

The primary purpose of analyzing the trajectory of the ERDAS-HYPACT and PROWESS-HYP ACT 
plumes was to determine how the model predictions compared with observations. Therefore, the part of the 
HYPACT plumes which were compared and vAnch are the focus of the discussion below are the lower plume 
below 925 meters and the upper plume above 1550 meters. Even though, the ERDAS-HYPACT simulation was 
run with the entire plume, we were not concerned with the transport of the center of the plume between 925 and 
1550 meters. 

3.233. HYPACT results 

The HYPACT predictions of the ERDAS and PROWESS plumes are shown in Figures 56 to 115. The 
figures show two-dimensional views of the plume from a horizontal map perspective and from a northward- 
facing vertical cross section. HYPACT output was produced eveiy ten minutes. 

HYPACT moved the lower cloud to the south vAxcrt it reached the coast in the vicinity of Cocoa Beach at 
approximately 30 minutes after the explosion. Figures 56 - 63 show the lowest part of the ERDAS cloud out 
over the ocean and extending upward and westward to the coastline where it curves north and then east as it 
extends upward. The cross-section shows the slight shift to the east of the particles in the lowest part of the 
plume. HYPACT begins spreading the lower cloud as shown by the diffusion of the lower plume particles 
(Figures 62 and 63). 

Figures 86 - 93 show the two PROWESS clouds (due to the initial source split) with the lower cloud to the 
southwest of the upper cloud. The lower cloud extends upward and westward from over the ocean to the 
coastline. The upper cloud extends eastward and upward from approximately 10 km offshore out to 
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approximately 20 km offshore. HYP ACT begins to also diffuse the lower PROWESS plume similar to the 
ERDAS plume. 

The observed radar plume between 1616 and 1705 UTC (Figures 17-23) show that the low-level cloud had 
reached the coast north of Cocoa Beach and had moved slightly more westward than the model runs had 
predicted. The orientation of the cloud was east-west in a shape similar to the model predictions. The upper 
cloud was also oriented in an east- west direction as it was transported eastward by the westerlies. 

Between 1710 and 1730 UTC, HYP ACT continued to diffuse both the ERDAS (Figures 64 - 69) and 
PROWESS (Figures 94 - 99) lower clouds as shown by the cross-sections which show considerable particle 
spread in the layers from the surface up to approximately 500 meters. HYP ACT continued to carry the lower 
cloud southward keeping most of the plume over the ocean extending barely on to the land. HYP ACT moved 
the ERDAS and PROWESS upper clouds rapidly to the east and causes very little spreading or diffusion within 
them. 

Figures 23 - 25, presenting the observed radar cloud from 1715 to 1735 UTC, show the lower cloud moved 
over the land further west than predicted by HYP ACT but kept the shape and orientation similar to that 
predicted by HYP ACT. The upper cloud exhibited an undulating shape as it moved south and east. 

During the period from 1740 to 1820 UTC, the ERDAS HYP ACT runs (Figures 70 - 75) and the 
PROWESS HYP ACT runs (Figures 100 to 105) moved the low level cloud to the south and spread it wider to 
the east and west. This trend continued for the remaining period of the simulation (Figures 76 - 85 and Figures 
110-115). The ERDAS low-level cloud extended slightly further west than the PROWESS low-level cloud. 
The upper cloud continued to move quickly to the east in both runs. 

From 1745 to 1833 UTC, the observed radar maps showed the movement of the two distinct clouds 
(Figures 26-31). The low-level cloud moved over Melbourne with part of the cloud extending eastward over 
the ocean. The bulk of the cloud was located over the land, to the west of where the model runs had predicted. 
The observed upper cloud matched the model predictions closely in its orientation from east-northeast to west- 
southwest and in its movement toward the east-southeast. Another interesting feature which was similar in the 
ERDAS HYP ACT upper plume (Figures 78 - 105) and the observed radar upper plume (Figures 24 - 36) was 
the wavelike, undulating shape. The wavelike structure formed in the direction of the strong westerly flow in the 
1000- to 1500-m layer. 

Over the final two hours of the simulation, from 1840 to 2030 UTC, the ERDAS (Figures 76 - 85) and 
PROWESS (Figures 106 - 1 15) HYP ACT runs continued moving the lower cloud south. Both runs diffused the 
lower clouds considerably with ERDAS HYP ACT causing more spread especially in the north-south direction. 

3.2.3.4. HYP ACT results summary 

One very interesting occurrence in the ERDAS HYP ACT simulation was the split that took place beginning 
at 1840 UTC between the upper and lower parts of the plume. The split occurred at height of approximately 
1 100 meters. During the final hour of the simulation from 1940 to 2030, the cross-section showed that the upper 
and lower part of the plumes looked distinctly different. The lower part was very diffuse and individual particles 
were visible. The upper part of the cloud was ribbon-like and showed little diffusion except for its elongation to 
the east as it was moved with the westerly flow. The importance of the decoupling of the plume was that the 
break occurred at the height of the inversion and it matched the two-cloud scenario which was observed visually 
and by radar. Both observed and modeled plumes had an undulating shape as discussed earlier. 

The discussion, provided in this section, indicates that the RAMS/HYPACT predictions of the lower 
cloud’s direction of movement, speed, and dispersion matched closely with observations. However, the onshore 
component of the winds affecting the observed plume appeared stronger than the RAMS modeled winds. As a 
result, the observed plume moved more to the west than the modeled plume indicated. The model’s prediction of 
the upper cloud’s direction and shape showed good agreement with the observed plume’s direction and shape. 

Comparing the PROWESS-HYP ACT plume to the ERDAS-HYPACT plume indicated that each model 
showed strengths and weaknesses when compared to the observations. 
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The source term used for the PROWESS-HYP ACT (split column source term) run was a better fit than the 
source term used for ERDAS-HYPACT (continuous column source term). This result is not sur^smg smce we 
decided to use the split column after viewing video and photography of the explosion cloud. The continuous 
column source term, however, showed a separation of the upper and lower clouds which ^gan to take place at 
approximately 1840 UTC at a height of about 1 100 meters. The separation was due to the strong shear above 
and below the strong inversion that existed and that was predicted by RAMS. 

The upper plume as predicted by ERDAS-HYPACT matched the observed upper plume slightly better than 
the upper pLne predicted by PROWESS-HYP ACT. PROWESS-HYPACT p am? 

than ERDAS-HYPACT and did not show the undulations that ERDAS-HYPACT predicted. e 
predicted winds, which control the plume movement in HYPACT, were predicted to be slightly stronger m 
PROWESS compared to ERDAS. 

At the surface, ERDAS-HYPACT and PROWESS-HYPACT both predicted approximately the same 
amount of spread of the lower plume in the east-west direction (Figures 79 and 107). The upper Pa« of the lower 
ERDAS-HWACT plume extended a little fiirther west than the PROWESS-HYPACT plume^ The obse^ed 
lower plume at 1853 UTC (Figure 33) was shifted slightly to the west of the modeled plumes and did not extend 
as far to the east over the water as predicted by the models. 
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Figure 56. ERDAS HYPACT, 1630 UTC, 17 Jan Figure 57. Same as Figure 56, vertical view from 
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Figure 58. ERDAS HYPACT, 1640 UTC, 17 Jan Figure 59. Same as Figure 58, vertical view from 
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Figure 60. ERDAS HYPACT, 1650 UTC, 17 Jan Figure 61. Same as Figure 60, vertical view from 
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Figure 62. ERDAS HYP ACT, 1700 UTC, 17 Jan Figure 63. Same as Figure 62, vertical view from 
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Figure 64. ERDAS HYP ACT, 1710 UTC, 17 Jan Figure 65. Same as Figure 64, vertical view from 
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Figure 66. ERDAS HYP ACT, 1720 UTC, 17 Jan Figure 67. Same as Figure 66, vertical view from 
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Figure 68. ERDAS HYPACT, 1730 UTC, 17 Jan Figure 69. Same as Figure 68, vertical view from 
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Figure 70. ERDAS HYPACT, 1740 UTC, 17 Jan Figure 71. Same as Figure 70, vertical view from 
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Figure 72. ERDAS HYPACT, 1800 UTC, 17 Jan Figure 73. Same as Figure 72, vertical view from 















7. = 157 fS m 


V =-200.43 krn 


1930 UTC 


Figure 80. ERDAS HYP ACT, 1930 UTC, 17 Jan Figure 81. Same as Figure 80, vertical view from 
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Figure 82. ERDAS HYPACT, 2000 UTC, 17 Jan Figure 83. Same as Figure 82, vertical view from 
I south. 
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Figure 84. ERDAS HYPACT, 2030 UTC, 17 Jan Figure 85. Same as Figure 84, vertical view from 
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Figure 86. PROWESS HYP ACT, 1630 UTC, 
17 Jan 1997. 
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Figure 88. PROWESS HYPACT, 1640 UTC, 
17 Jan 1997. 
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Figure 90. PROWESS HYPACT, 1650 UTC, 
17 Jan 1997. 


y =-280.43 km 


1630 UTC 


Figure 87. Same as Figure 86, vertical view from 
south. 
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Figure 89. Same as Figure 88, vertical view from 
south. 
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Figure 91. Same as Figure 90, vertical view from 
south. 
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Figm-e 92 PROWESS HYPACT, 1700 UTC, Figure 93. S«ne us Figure 92, vertical view from 
17 Ja» 1997. 



Figure 94. PROWESS HYPACT, 1710 UTC, 
17 Jan 1997. 
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Figure 95. Same as Figure 94, vertical view from 
south. 
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Figure 96. PROWESS HYPACT, 1720 UTC, 
17 Jan 1997. 
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Figure 97. Same as Figure 96, vertical view from 
south. 
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Figure 98. PROWESS HYP ACT, 1730 UTC, 
17 Jan 1997. 
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Figure 100. PROWESS HYP ACT, 1740 UTC, 
17 Jan 1997. 
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Figure 102. PROWESS HYP ACT, 1800 UTC, 
17 Jan 1997. 
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Figure 99. Same as Figure 98, vertical view from 
south. 
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Figure 101. Same as Figure 100, vertical view from 
south. 
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Figure 103. Same as Figure 102, vertical view from 
south. 
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Figure 104. PROWESS HYP ACT, 1820 UTC, 
17 Jan 1997. 
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Figure 106. PROWESS HYPACT, 1840 UTC, 
17 Jan 1997. 
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Figure 108. PROWESS HYPACT, 1900 UTC, 
17 Jan 1997. 
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Figure 105. Same as Figure 104, vertical view from 
south. 
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Figure 107. Same as Figure 106, vertical view from 
south. 
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Figure 109. Same as Figure 108, vertical view from 
south. 
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Figure 110. PROWESS HYP ACT, 1930, UTC, Figure 111. Same as Figure 110, vertical view from 
17 Jan 1997. south. 
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Figure 112. PROWESS HYP ACT, 2000 UTC, Figure 113. Same as Figure 112, vertical view from 
17 Jan 97. south. 



Figure 114. PROWESS HYP ACT, 2030 UTC, Figure 1 15. Same as Figure 1 14, vertical view from 
17 Jan 97. south. 



3.2.4. ERDAS HYP ACT trajectory analysis 


We conducted an analysis to determine how well the model predicted the track of the plume when 
compared to observations. We did not compare observed versus predicted concentrations because no 
concentration measurements were available for this case. However, we used concentration predictions from 
ERDAS HYP ACT to help identify the model’s prediction of plume location. HYP ACT was configured to 
output concentrations every 10 minutes and in layers with a thickness of 75 meters, HYP ACT produced 
concentrations on a 3-dimensional grid at 10-minute intervals from 1630 UTC until 2030 UTC. The location of 
the observed plume was determined from the WSR-88D analysis. The location of the predicted plume was 
determined by finding the maximum concentrations predicted by the ERDAS HYP ACT model for each 10- 
minute period in each vertical layer of interest. The latitude, longitude, and height of this maximum 
concentration were then noted. By tracking this maximum at the 10-minute intervals we were able to determine 
the trajectory of the HYPACT-predicted plume. 

We had determined from the visual and radar observations of the initial explosion plume that its 
dimensions were not precisely defined. However, the plume extended vertically from the ground up to a height 
of approximately 900 meters. Because we did not know the height of maximum concentrations, we chose a 
layer that was representative of the lower plume. We chose a layer centered at 412.5 meters since this layer was 
near the midpoint of the lower plume. 

The graph in Figure 116 shows the location of the maximum concentration of the layer centered on 412.5 
m. This layer is near the center of the ground-based layer which contained the lower plume resulting from the 
explosion. The plume originated at Launch Complex 17 and then moved south over the following 4 hours. The 
observed and predicted plume locations closely agreed for the first two hours except that the observed plume 
was slightly west of the predicted plume. During the last two hours of the model simulation, the graph gives the 
appearance of a widely spreading plume. This spread was not due to the plume shifting from the east to the west 
but due to the shift in the location of the predicted maximum concentration within the plume from one 10- 
minute period to the next. 
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Figure 1 16. Data of observed vs. predicted plume location for lower plume resulting from Delta II explosion at 
1628 UTC on 17 January 1997. The observed plume location was determined by WSR-88D radar. The 
predicted plume location was determined by finding the maximum predicted ERDAS HYP ACT concentration 
in the 75-m layer centered at a height of 412 m. The numbered circles represent the time of the observed and 
predicted plumes with the following time tags: 1: 1600-1700 UTC, 2: 1700-1800 UTC, 3: 1800-1900 UTC, 4: 
1900-2000 UTC, 5: 2000-2100 UTC. 
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4. Conclusions and Recommendations 

The Drimarv goal of this study was to conduct a case study of the dispersing cloud and the tnodels used to 

pi J= obse^aaops. The „«h»rologiea, eo„«o„s on hhs day were singly 
influenced by synoptic rather than local forcing. 

The conelnalons of this study ean be eategoriaed aeeorfing to the plnn« obsetvation technique and 
according to the models used in the analysis. The findings of this smdy are. 

WSR-88D radar as a plume observation tool 

. T t,. WS^R.R8D is a vood tool for prov i Hina nlume tracks from rocket explosion pte The rad^ 
provided excellent data on the location and tock of the resultmg potentially toxic cloud. The data 
was extremely useful for model verification since no ongoing program is m place to measure 
plume track or concentrations. Bud Parks of ACTA, Inc. is conducting a study to ^ 

nominal and abort launch clouds and has captured ^ data Eastern Range abort 

clouds and 36 of 46 (78%) nominal launch clouds (Parks and Evans 1998). 

• The WSR-88D does not provide concentration date. The only data obtained by the radar is the 
reflectivity value measured in dBZ. While this measurement gives an estimate on &e relative 
delity of material (smoke particles, water, dust, or other particulate matter), a methodology is 
needed to convert dBZ to concentrations of hydrochlonc acid (HCl), mtrogen tettoxide (N2O4), o 
other materials of interest. One of Range Safety’s main concerns is determining the exposure limit 
(concentration over a specified time) of certain toxic materials. 

A dark orange cloud at the very top of the large lower cloud w^ initially visible. The tok oi^ge 
cloud most ILkely contained some amount of nitrogen tetroxide. Bemuse it 

of the cloud, it is unsure how much, if any, of the N2O4 mixed witto Ae cloud and made it to the 
surface. Our analysis was not able to determine the concentration of N2O4 in the explosion cloud. 

. nlnme height data for this case was not v ery accurate. The radar appeared to accurately 

track the clouds’ trajectories in the x-y dimension. However7*e vertical me^urements appe^d 
to be inaccurate for two reasons. The first reason was that for long distances from &e such 

as the 35+ kilometers distance from Melbourne to Cape Canaveral, the radar beam widened 
enough to introduce inaccuracies in the vertical plume height measurements. 

inaccurate vertical measurements was because of the strong mvers.on causmg the radar beam to 
bend and bring about measurement inaccuracies. 

RAMS model 

• The vertical wind profile predicti o ns in this case show agreement with observa^on s. Both E^ AS 
and PROWESS configur^ions of RAMS produced wind flow measurements that niatc^d closely 
with rawinsonde and profiler measurements and seemed to provide adequate mput to HYP ACT. 

. p Atvi. nnHcr-nredicteH on.bore flow at the level o f the Delta II .cloud. Bodi 

predicted onstore flow in the 600- to 900-meter layer m the area south ^P® 

Lwever the observed winds had more of a northeast, onshore component th^ the RAMS- 
predicted winds. The movement of the actual explosion cloud tracked onshore north of the model- 
Jredicted cloud. The prototype ERDAS RAMS predicted wmds with more of an onshore 

component than PROWESS RAMS. 

. RAMS . ccnratelv Predicted the strength a nd height of inversion for this case. The well-defined 
inversion that iJs measured by rawiLonde and based at approximately 800 meters above the 
ground was predicted by RAMS to be based at approximately 750 meters above the ground. The 


59 


inversion, as determined by the vertical temperature profile, had a significant influence on the 
explosion cloud. 

R££DM source term 

• Characterizing the source term of unique explosions is difficult. If a rocket explosion occurs, the 
circumstances will be different each time it happens. For example: What was the flight time? How 
much fuel was consumed? What were the height, location and distribution of the explosion 
products? Were the hazardous and toxic materials separated or mixed within the cloud? Did the 
second stage ascend and then explode as with the Delta II? Did the solid rocket motors explode 
immediately or did they follow an errant path before they exploded? All of these questions make it 
difficult to develop a model that will accurately assess and characterize the source term. We were 
able to use information obtained after-the-fact from radar and video to characterize the source term 
but in a real-time scenario only estimates of the source term characterization can be made. 

• Splitting the source into two sources for HYP ACT model was a reasonable approximation. In the 
post-explosion mode of this analysis, we split the source into two sources for the HYP ACT model 
based on observations. This methodology proved better for this case as opposed to using the single 
column source term that was generated by REEDM. In real-time, REEDM currently does not split 
the source term. 

HYP ACT model 

• Plume came onshore further to the north and earlier than predicted. HYP ACT moved the large 
lower cloud resulting from the Delta II explosion onshore at a point that was approximately 12 
kilometers south of where it actually came onshore. HYP ACT predicted the plume would come 
across the coastline in the Satellite Beach/Indian Harbor Beach area when it actually crossed the 
coastline in the Cocoa Beach area. The observed winds at the level moving the cloud at 600-900 m 
were northerly prior to the explosion but then shifted to northeasterly during the 1-2 hours 
following the explosion. RAMS predicted the shift of these winds from northerly to northeasterly 
but did not predict the shift as quickly as observed. RAMS supplies wind and meteorological data 
to HYP ACT. Because of RAMS* gradual response to shift the winds, HYP ACT missed the 
location and timing of the plume impact on the coastline. 

• Trajectory, diffusion and timing of HYP ACT plumes showed similarities to observed plume. 
Except for the problem mentioned above, the trajectory, diffusion, and the timing of the HYP ACT 
plumes were similar to the plumes observed by radar. One favorable result was noted in the spread 
and diffusion of the lower cloud as it moved south. The cloud spread in the crosswind direction at 
a rate and in distance similar to observed. 

• Range Safetv*s REEDM predicted the movement of plume to the south (176 degrees). REEDM 
using the 1613 UTC rawinsonde from Cape Canaveral moved the plume to the south and kept it 
offshore until it reached the Melbourne Beach area. REEDM did not account for the winds with an 
easterly component that existed at a height of 700-800 meters in the area over the ocean to the 
south of Cape Canaveral. 

Reconunendations 

• Develop methodology to correlate concentrations with rztdar reflectivity measurements. The WSR- 
88D proved to be a valuable tool in tracking nominal and abort rocket plumes. However, the radar 
provides no information on the concentrations within the clouds. What is needed is measurement 
of concentrations within the plumes using a sample collection method or another remote sensing 
technique such as lidar. This data could then be correlated with radar measurements of reflectivity 
indBZ. 

• Improvements are needed in HYP ACT plume dynamics algorithms. HYP ACT currently treats 
plumes as non-buoyant, non-depositing entities. We recommend that future enhancements should 
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be made to HYP ACT to improve its ability to handle buoyant plumes and particle deposition. 
These improvements would allow HYP ACT to model rocket exhaust plumes better than the 
current version of HYP ACT. 


Conduct other studies of roc ket explosion plumes. Since the explosion of the Delto H two other 
rockets have exploded after launch from Cape Canaveral-Titan IV on 12 August 1998 ^d Delta 
III on 26 August 1998. In both cases the explosion clouds were tracked by WSR-88U radar. 
Detailed studies should be conducted to verify mesoscale models, diflusion models, and radar 


tracking techniques. 
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APPENDIX A 
ERDAS RAMSIN file 
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$MODEL_GRIDS 

I Simulation title (64 chars) 

EXPNME = 'CCAFS/KSC -Version 3a’ , 


lOTYPE 

0. 



runty PE = 

•INITIAL’ , 


TIMEUNIT = 

•h' , 



TIMMAX= 

12. , 



ngrids = 

3, 



NNXP 

38, 

34, 

37, 

NNYP = 

36, 

38, 

37, 

nnzp 

25, 

25, 

25, 

NNZG 

11, 

11, 

11, 

NXTNEST = 

1, 

1, 

2, 


» - l=diagnostic , 0=normal run 

type of run; INITIAL, HISTORY, or MAKEVFILE 

1 ’h'r’m',’s' - Time units of 

i TIMMAX, TIMSTR. VTIME 

! Final time of simulation 

! Number of grids to run 

I 

! Number of x gridpoints 

! Number of y gridpoints 

! Number of z gridpoints 

! Number of soil gridpoints 

1 Grid number which is the next 
i coarser grid 


! Coarse grid specifications 


IHTRAN = 
DELTAX = 
DELTAY = 

1, 

60000. , 
60000. , 


! 0-Cartesian, 1-Polar stereo 
! X and Y grid spacing 

DELTAZ = 

DZRAT 

DZMAX 

75. , 
1.25, 
1000. , 


! Z grid spacing (set to 0 . to use ZZ) 

» vertical grid stretch ratio 
1 maximum delta Z for vertical stretch 

ZZ=0 . , 25 - , 75 . , 150 . 
3250., 4000., 5000 

, 250 . , 500 . , 
6000., 

750., 1000., 1500., 2000., 2500., 

I Vertical levels if DELTAZ = 0 

DTLONG = 
NRATIO = 

90. , 
3, 


! Coarse grid long timestep 
1 Small timestep ratio 

IMONTHl = 
IDATEl = 
lYEARl = 
STRTIM = 

01, 

17, I 

97, 

12.0. 

! Month 
Day 
Year 

! GMT of model TIME = 0. 




! Nest ratios between this grid 
i and the next coarser grid. 

NSTRATX = 
NSTRATY = 
NNDTRAT = 

1, 4, 5, 
1, 4, 5, 
1, 2, 2, 


1 x-direction 
1 y-direction 
! time 

NESTZ = 3, 
NSTRATZ=3, 

3, 2, 2, 2, 2, 

2,2, 2, 2,1, i 

1 contort coarser grids if negative 

POLELAT 
POLELON = 

= 30.0, 

-83.0, 


! Latitude of pole point 
! Longitude of pole point 

CENTLAT= 

CENTL0N= 

30.0, 

-83.0, 


1 Grid point on the next coarser 
1 nest where the lower southwest 
' corner of this nest will start. 

! IF NINEST or NJNEST = 0, use CENTLAT/LON 
! coarser grid 

NINEST = 
NJNEST = 
NKNEST = 

1, 17, 
1, 10, 
1, 1. 

20, 

19, 

1, 

i i -point 
t j -point 
I k-point 

NNSTTOP = 
NNSTBOT = 

1, 1, 

1, 1. 

0, 

1, 

! Flag (0-no or 1-yes) if this 
! nest goes the top or bottom of the 
! coarsest nest. 

GRIDU 

gridv 

0. , 0. , 

0., 0., 

0. , 
0. , 

* u- component for moving grids 
! V- component for moving grids 
! (still not working!) 
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$END 


$MODEL_FILE_INFO 


I History file input 

TIMSTR= 0., ! time of history start (see TIMEUNIT) 

HFILIN = 'ci.hOh', ! input history file name 

! History/analysis file output 

IOUTPUT= 1, ! 0-no files, 1 - save/dispose 

HFILOUT= ’ hist/f .h ’ , ! history file prefix 

AFILOUT= ’ f . a ’ , ! analysis file prefix 

HFUNITS=’h', AFUNITS='d', I history/anal file units (M,m, H, h, S, s) 

FRQHIS =43200., FRQANL =3600., ! history/anal file frequency 

I Variable initialization input 

INITIAL =2, 1 Initial fields - l=horiz .homogeneous , 2=variable 

VTIME = 0., 6., 12., 18., 24., I model times for varfiles (See TIMEUNIT) 

VARFIL(l) = ’isan/ivl7- jan-97-12' , VTIME(1)=0. , 

VARFIL(2) = *isan/ivl7- jan-97-18’ , VTIME (2) =6. , 

VARFILO) = ’ isan/ivl7- jan-97-24 ’ , VTIME (3) =12 . , 

VARFIL(4) = ' isan/ivl7- jan-97-30' , VTIME (4) =18 . , 

VARFIL(5) = 'isan/ivl7- jan-97-36' , VTIME (5) =24. , 

IVWIND = 0, ] initial winds ( only 0 works) 

NUDLAT =5, ] number of points in the lateral bnd region 

TNUDLAT = 900., I nudging time scale (s) at lateral boundary 

TNUDCENT = 00., ! nudging time scale (s) in center of domain 

TNUDTOP = 00., ! nudging time scale is) at top of domain 

ZNUDTOP = 15000., ! nudging at top of domain above this height (m) 

! Printed output controls 

FRQPRT = 21600., ! Printout frequency 

FRQIPR = 99930., ! Integral print frequency 

FRQIST = 99900., ! Integral store frequency 

ISTPFL = 1, i Timestep message frequency flag 

INITFLD = 1, ! Initial field print flag 0=no print, l=print 

INPRTFL = 1, ! Namelist print flag 0=no print, l=print 

! Input topography variables 

SFCFILES = • /ul/met/sf c/sf cb _jan', I File path and prefix for surface 

characteristic 

ITOPTFLG =1,1,1, I 2 - Fill data in «ruser*' 

IPCTLFLG =1,1,1, ! 1 - Interp data from latlon dataset 

ISSTFLG =1,0,0, ! 0 - Interpolate from coarser grid 

IVEGTFLG = 1,1,1, 


! The following only apply for IxxxxFLG=l 

; 

ITOPTFN = ’ /usr/geode/data/rams/topolOm/H ‘ , 

' /usr/geode/data/rams/topo30s/U' , ! Input topography file name 

' /usr/geode/data/rams/topo30s/U’ , ! Input topography file name 

IPCTLFN = ' /usr/geode/data/rams/pctllOm/L' , 

’ /usr/geode/data/rams/pctllOm/L’ , 

• /usr/geode/data/rams/pctllOm/L* , 

ISSTFN = ' /usr/geode/data/rams/sst/SJUL ' , 

IVEGTFN = ’ /usr/geode/data/rams/ndvi/V , 

’ /usr/geode/data/rams/ndvi/V , 

‘ /usr/geode/data/rams/ndvi/V , 


SILAVWT = 0.0, 0., 0., ! Weighting of topo silhouette averaging 
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toptwvl 

= 4.0, 

4 

PCTLWVL 

= 2.0, 

2 

SSTWVL 

= 2.0, 

2 


! Topo wavelength cutoff in filter 
1 Land pet wavelength cutoff in filter 
! Land pet wavelength cutoff in filter 


$END 


$MODEL_OPTIONS 
NADDSC = 0,0,0, 


NTOPSMTH 
I Z FLAT = 


0, ! Number of passes for 

o[ \ Width of flat margin 


topography smoother 

around domain (in grid points) 


1 Numerical schemes 


ITMDIFF 

NONHYD 

SSPCT 

IMPL 

ICNTEQ 

WTKD 

ICORFLG 


3, 
1 , 
0 . , 
1 , 
2 , 
0 . , 

1 , 


l=forward, 2=leapfrog, 3=hybrid 
nonhydrostatic=l, hydrostatic=0 
Sound speed fraction for the nonhydros ta tic 
implicit flag for acoustic model - 

Hyd - continuity equation - 1 -incomp, 2 -anelastic 
Klemp/Durran - current times tep weighting 






. p>n t" 




IBCTOP = 


IBND 

JBND 

CPHAS 

LSFLG 


NFPT 

DISTIM = 

TIMSCL = 
KSPIN 
KMSPIN = 

IPRSPLT = 
lADVL 
lADVF 
IPGRAD = 

FILT4 


FXLONG = 
FYLONG = 


top boundary condition ^ 

0- wall on top(nonhyd) 2 -prognostic sfc 

1 - Klemp-Durran (nonhyd/hyd) 3 -material surface (hyd) 


1 , 

1 . 

20 . , 

0 , 


0 , 

0 . , 

0 . . 
33, 
10 , 

0 , 

2 , 

2 , 

1 , 

000 . , 


0 . 0000 , 

0 . 0000 , 


3 -Orlanski 


! Lateral boundary condition flags 
! 1 - Klemp/Wilhelmson, 2-Klemp/Lilly, 

! 4 -cyclic 

! Phase speed if IBND or JBND - 1 

! Large-scale gradient flag for variables other than 
! normal velocity; 

! 0 = zero gradient inflow and outflow 
! 1 = zero gradient inflow, radiative b.c. outflow 
! 2 = constant inflow, radiative b.c. outflow 
! 3 = constant inflow and outflow 

Rayleigh friction - number of points from the top 
I - dissipation time scale 


Initial wind spin-up 


time scale 

- below this level (velocity) 

- below this level (momentum) 

precipitation time- split scheme - 
Order of advection - Leapfrog - 2 or 4 
Order of advection - Forward - 2 or 6 
pressure gradient scheme for topo.: 

Fourth order filter coefficient (0 - off) 

{ > 0 times tep at which 2 delta waves are 
totally removed) 

Long filter coefficient 0=off 
Long filter coefficient 0=off 


! Radiation parameters 


ISWRTYP = 2, 
ILWRTYP = 2, 
RADFRQ = 1200. 
LONRAD = 1 / 


Shortwave radiation type ^^Mahrer/Pielke, l_Chen 

Longwave radiation type 0-none, 2-Mahrer/Pielke, 1 Che 
I Frequency of radiation tendency update in seconds 
Longitudinal variation of shortwave (0-no, 1-yes) 


I Cumulus parameterization parameters 


NNQPARM = 0, 0, 

CONFRQ = 1200., 
WCLDBS = .001, 


I ! convective parameterization flag 

! 0-off, 1-on 

'reauency of conv param. updates in seconds 
vertical motion needed at cloud base for convection 


! Surface layer and soil parameterization 


ISFCL = 


NVGCON = 


1, ! surface layer/soil/veg model 

1 0- specified surface layer gradients 
! 1 - soil /vegetation model 

I, ! Vegetation type 

1 1 - - Crop/mixed farming 


Short grass 
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! 3 - - Evergreen needleleaf tree 4 - - Deciduous needleleaf tree 

! 5 - - Deciduous broadleaf tree 6 - - Evergreen broadleaf tree 

! 7 -- Tall grass 8 -- Desert 

! 9 - ■ Tundra 10 - - Irrigated crop 

ill -- Semi-desert 12 -- Ice cap/glacier 

113 -- Bog or marsh 14 -- Inland water 

! 15 - - Ocean 16 - - Evergreen shrub 

!17 -- Deciduous shrub 18 -- Mixed woodland 

TSEASN = 298., ! Average seasonal temp 

TVGOFF = 0-, ! Initial veg temp offset 

VWTRCON = .000, ! Initial veg water storage (not working) 

UBMIN = 0.25, ! Minimum U value to use in computing 

PCTLCON =1.0, ! constant land percentage if for all domain 

NSLCON =6, i constant soil type if for all domain 

! l=sand 2=loamy sand 3=sandy loam 

! 4=silt loam 5=loam 6=sandy clay loam 

! 7=silty clay loam 8=clay loam 9=sandy clay 

! 10=silty clay ll=clay 12=peat 

ZROUGH =0.25, ! constant roughness if for all domain 

ALBEDO = 0.2, I constant albedo when not running soil model 

SEATMP = 280., I constant water surface temperature 

DTHCON = -10., 1 constant surface layer temp gradient for no soil 

DRTCON = .000, ! constant surface layer moist gradient for no soil 

SOILDZ =0., ! soil model grid spacing 

SLZ = - . 50, - .40, - . 30 , - . 25 , - . 20 , - . 16 , - . 12 , ' . 09, - . 06 , - . 03 , 0.,! soil grid levels 

SLMSTR = 0.6, 0.6, 0.5, 0.5, 0.5, ! initial soil moisture 

0.5, 0.5, 0.4, 0.4, 0.4, 0.4, 

STGOFF= 5., 5., 5., 5., 3.5, 2., .5, -1., -1.5, -1.8, -2., 

I Initial soil temperature offset from lowest 
I atmospheric level 


! Eddy diffusion coefficient parameters 

IDIFFK =1,1,1, » K flag: 

I 1 - Horizontal deformation/ Vertical Mellor- Yamada 

I 2 - Anisotropic deformation (horiz and vert differ) 

! 3 - Isotropic deformation (horiz and vert same) 

I 4 - Deardorff TKE (horiz and vert same) 

CSX = .2, ,2, ,2, ! momentum M-Y k 

CSZ = 20. ,20. ,20., ! scalar M-Y K 

XKHKM = 3.,3.,3., ! Ratio of horizontal K_h to K_m for TKE or deformation 

2KHKM = 3.,3.,3., ! Ratio of vertical K__h to K_m for TKE or deformation 

AKMIN= .3,. 3,. 3, ! Ratio of minimum horizontal eddy viscosity coefficient 

I to typical value from deformation K 


Microphysics 


NLEVEL 

= 

1, 

1, 1 

! moisture complexity level 

I CLOUD 

= 

0, 

0, 0, 



IRAIN 

= 

0, 

0, 0, 

i Microphysics 

flags 

IPRIS 

= 

0, 

0, 0, 



ISNOW 

= 

0, 

0, 0, 

5 where x= R - 

rain 

lAGGR 

= 

0, 

0, 0, 

! P - 

pristine crystals 

IGRAUP 

= 

0, 

0, 0, 

S S - 

snow 

IHAIL 

= 

0, 

0, 0, 

! A 

aggregates 

CPARM 


0. . 

0. ,0. , 



RPARM 

= 

0. , 

0. ,0. . 

1 G - 

graupel 

PPARM 

= 

0. , 

0. ,0. , 

! H - 

hail 

SPARM 

= 

0. , 

0 ., 0 ., 

! NIxCNFL = 0 

no species 

APARM 

= 

0. , 

0. ,0. , 

1 1 

diagnostic concen . 

GPARM 

= 

0. , 

0. ,0. , 

! 2 

specified mean diameter 

HPARM 

= 

0. , 

0. .0. , 

! 3 

specified y- intercept 





! 4 

specified concentration 





! 5 

prognostic concen. 

AMIO 

- 

0.0000001, ! 

minimum crystal mass (gm) 


CON 


300 

.0 , ! 

CCN concentration (number/cc) 

THOMO 

= 

233. 

0, ! 

homogeneous nucleation temperature 


$END 

$MODEL_SOUND 

! Sounding specification 
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IPSFLG=0, 


! Flags for how sounding is specified 
! specifies what is in PS array 

! 0 - pressure (inb) , 1 ■ heights (m) , PS(l)=sfc press (n\b) 

ITSFLG=2, I specifies what is in TS array 

! 0 - temp(C), 1 - temp(K), 2 - pot. temp(K) 

IRTSFLG=3, I specifies what is in RTS array , , 

! 0 ' dew pnt. (C) , 1 - dew pnt. (K) , 2 - mix rat(g/K.g) 

! 3 - relative humidity on %, 4 - dew pnt depression (K) 

IUSFLG=0, ! specifies what is in US and VS arrays 

! 0 ' u,v component (m/s ) , 1 ■ umoms- direction, vmoms* speed 

lUSRC = 0, ! source of wind profile: 

! 0 - umoms, vmoms, valid at sounding levels (PS) 

! -1 - usndg, vsndg, valid at model levels (Z) 


PS=1012 . 00. 

, 1000.00 

, 977.00 

, 949.00 

, 936.00 

, 850.00 

, 848.00 

797.00, 

700.00 , 

662.00 , 

595.00 , 

579.00 , 

565.00 , 

560.00 , 

500.00, 

498.00 , 

461.00 , 

443.00 , 

406.00 , 

400.00 , 

387.00 , 

383.00, 

328.00 , 

325.00 , 

316.00 , 

309.00 , 

307.00 , 

300.00 , 

285.00, 

250.00 , 

200.00 , 

150.00 , 

100.00 , 

80.00 , 

10.0, 


TS= 

288.18, 

288.36 , 

288.88 

293.11 , 

294.27 , 

298.93 , 

298.92 

304 .48, 

309.35 , 

310.27 

313.50 , 

315.02 , 

317.23 , 

317 . 80 


318.75, 

318,88 , 

320.88 

322.03 , 

326.27 , 

327.02 , 

328.41 


328.73, 

330,56 , 

331.02 

333.69 , 

335.55 , 

335.89 , 

336.99 


340.82, 

343.42 , 

347.35 

364.39 , 

379.82 , 

398.66 , 

480. , 

RTS^ 

= 90.7, 

95.5, 

98.0, 

96.8, 

96.1, 

100.0, 

100.0, 

94.1, 

90.6, 

99.3, 

94.1, 

75.8, 

49.3, 

38.4, 

55.8, 

55.7, 

85.1, 

80.6, 


78.0, 

09.7, 

71.1, 

02.5, 

57.6, 

06,1, 

52.3, 

15.8, 

74.9, 

03.6, 

74.0, 

02.1, 

39.6, 

01.5, 

22.7, 

01.5, 

16.0, 


US=4 ., 4 . ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4., 

4 . " 4 . ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4. ,4., 

VS=0 0 ,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0w0.,0.,0.,0., 

0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ,, 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 

US= -10., -10., -11., -12., -12., -5., 5., 

6.5, 8.7,11.0,14.0,16.0,17.0,18.0, 

21.9, 22.0,24.0,26.0,30.0,30.0,32.0, 

33.0, 34.0,35.0,36.0,37.0,39.0,40.6, 

42.0, 50.7,46.6,35.0,30.0,30.0,30., 

VS= 2.1, 2.1, 2 . 5 , 3. 0,3. 5, 4. 0,4.0, 

3.5, 3.1, 3 . 3 , 3 . 5, 3. 7, 4. 2, 4. 5, 

6.0. 6.0, 7.0,8.0,9.0,10.0,11.0, 

! 12.0. 13.0,14.0,14.1,14.5,15.0,15.0, 

1 14.8, 8 . 0 , 8 . 0 , 8 . 0 , 8 . 0,8. 0,8.0, 

USNDG =0,, ! Wind components if IUSRC= -1 

VSNDG = 0 . , ! 

KMEANl =0, ! lower model level for calculation of (umean, vmean) 

KMEAN2 =0, ! upper model level for calculation of ( umean , vmean ) 

UMEAN =0.0, 1 u- component for Galilean transformation 

VMEAN =0.0, I v- component for Galilean transformation 

$END 


$MODEL_PRINT ^ 

I Specifies the fields to be printed during the simulation 

NPLT = 6 , ! Number of fields to be printed at each time 

! for various cross-sections (limit of 50) 
IPLFLD = 'UP' , 'VP' , 'WP’ , ’THETA' , ’ TKE ’ , 'TVEG' , 

! Field names - see table below 

! PLFMT( 6 ) = ’0PF7.2', ! Format specification if default is unacceptable 
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IXSCTN = 3 , 3 , 3 , 3 , 3 , 3 , 3 , 3 , 3 , 3 , 3 , 3 , 

! Cross ' section type (1— XZ, 2-YZ/ i-XY] 

ISBVAL = 2, 2,2,2, 2, 2, 2. 2, 2, 2, 2, 

I Grid-point slab value for third direction 

< The following variables can also be set in the $PRNT namelist: lAA, 

! lAB, JOA, JOB, NAAVG, NOAVG, PLTIT, PLCONLO, PLCONHI , and PLCONIN. 


$END 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 


'UP' 

• VP' 

'WP' 

' PP' 

'THP' 

’ THETA • 

»THVP’ 

'TV 

' RV 

'THIL' - 

• TVP • 

' RELHUM ' 
' FTHRD • - 
'MICRO' - 
■ ZO ' 

' USTARL ’ 
• TSTARW ' 
'UW' 

’ WFZ • 
'QFZ' 

' RSHORT' 


- UP(M/S) 

'RC 

- RC(G/KG) 

•PCPT' 

- TOTPRE 

- VP(M/S) 

'RR' 

- RR(G/KG) 

'TKE' 

- TKE 

- WP(CM/S) 

'RP' 

- RP(G/KG) 

'HSCL* 

- HL(M) 

- PRS(MB) 

- THP{K) 

'RA' 

■ RA{G/KG) 

'VSCL* - 

VL(M) 

- THETA (K) 

• RL' 

■ RL(G/KG) 

'TG' 

TG (K) 

- THV* (K) 

'RI* 

• RKG/KG) 

'SLM* 

SLM (PCT) 

- TV(K) 

• RCOND ' 

- RD(G/KG) 

■ CONPR ' - 

CON RATE 

- RT<G/KG) 

'CP' 

- NPRIS 

' CONP • - 

CON PCP 

- RV(G/KG) 


- RT' {G/KG) 

• CONH ' - 

' CONM • - 

CON HEAT 
CON MOIS 


Theta-il (K) 'TEMP' - temperature (K) 

1V« (K) 'THV - Theta-v (K) 

-relative humidity (%) 'SPEED'- wind speed (m/s) 

radiative flux convergence (] 

GAS PRC 

ZO (M) 'ZI' - ZI (M) 

-USTARL(M/S) 'USTARW -USTARW(M/S) 

-TSTARW(K) ' RSTARL' -RSTARL (G/G) 

UW (M*M/S*S) 

WFZ (M*M/S*S) 

QFZ (G*M/G*S) 

- RSHORT 


' ZMAT ’ - 

' TSTARL ' 
’ RSTARW ' 
• VW' 
•TFZ' 

’ RLONG ' - 


ZMAT (M) 
TSTARL (K) 
RSTARW (G/G) 
VW {M*M/S*S) 
TFZ (K*M/S) 
RLONG 


$ISAN_CONTROL 
MSTAGE = 1 , 


NATIMES = 1, 

lAHOUR = 00, 
I ABATE = 18, 
lAMONTH = 7, 
lAYEAR =91, 


1 , 


! Main switches for 

! pressure, isentropic, "varfile" processing 
Number of times on which to perform analysis 


06, 

12, 

18, 

24, 

18, 

19. 

19, 

20, 

7, 

7, 

7, 

7, 

91, 

91, 

91, 

91, 


! Hours to analyze 
I Dates 


i Months 
! Years 


! NCAR archived data file names 
lAPR = 'grid.O', 5 Input pressure level dataset 

lARAWI = ' tempdb ' , ^ Archived rawindsonde file name 

lASRFCE = 'not -used' , ! Archived surface obs file name 

! File names and dispose flags 


IFNPRS = ’ ip' , 
IOFLGP= 1, 

IFNISN = ' ii ' , 
IOFLGI= 1, 

IFNSIG = * is* , 
IOFLGS= 1, 

IFNVAR = ' iV , 
IOFLGV= 1, 


I Pressure file name prefix 

1 Dispose flag: 0 = no write, 1 = write 

! Isentropic file name prefix 

! Dispose flag: 0 = no write, 1 = write 


! sigma- z file name prefix 

! Dispose flag: 0 


no write, 1 = write 


"Variable initialization" file name prefix 

1 Dispose flag: 0 = no write, 1 = write 


$END 

$ISAN_PRESSURE 


Pressure grid information: 
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NPRX=36, 

NPRY=33, 

NPRZ=10, 


! number of grid points in x (Ion) direction 
! number of grid points in y (lat) direction 
! number of pressure levels 


PDATA= 'SPC, 

WPLON = -140.0, 

SPLAT = 20.0, 


! 'NMC or ’ECMWF' if from NCAR archives. 

I West longitude bound of pressure data access 
! South latitude bound 

! Latitude and longitude bounds are north and east 
! if positive, south and west if negative. 


SPCNPRX = 2.5, 
SPCNPRY = 1.25, 


! East -west grid spacing of pressure data 
1 North- south spacing 


! Pressure levels (mb) in input dataset to access 

LEVPR = 1000,850,700,500,400,300,250,200,150,100, 


$END 


Isentropic and sigma* z processing 


$ I SAN_ISENTROPIC 


! Specify isentropic levels 


NISN 

LEVTH 


32, ! number of isentropic levels 

282 284 286 288,290,292,294,296,298,300,303,306,309,312, 

3 ? 5 ; 318 : 32 ^ 324 : 327 ; 330 ; 335 , 340 , 345 , 350 , 3 ^ 

440,460, 


NIGRIDS = 

TOPSIGZ 

HYBBOT 
HYBTOP 

SFCINF 

SIGZWT 


= 20000 ., 

= 4000. , 

= 6000, . 

= 1000 . , 

= 1 . , 


Analyzed grid information: 

1 number of RAMS grids to analyze 
! sigma* z coordinates to about this height 


NFEEDVAR = 


- Bottom (m) of blended sigma -z/isentropic layer in varfiles 
! Top (m) of blended sigma- z/isentropic layer 

! vertical influence of surface observation analysis 

! Weight for sigma- z data in varfile: wvpmoT 

! 0.= no sigz data, l.=full weight from surface to HYBBOT 

1 1 = feed back nested grid varfile info, 0 = don't 


Observation number limits : 


MAXSTA=1S0, 

MAXSFC=1000, 

lARCSND = 0 , 

lARCSFC = 0 , 

ISPCSND = 1, 

ISPCSFC = 1/ 

NONLYS = 0 , 

IDONLYS =’76458 ’ , 

NOTSTA = 0 , 

NOTID =’r76458' 


STASEP = .0001, 


ISTAPLT = 0, 


I maximum number of rawindsondes (archived + special) 
! maximum number of surface observations 


! Input archived soundings ? (0-no, 1-yes) 

! Input archived surface obs ? 

! Input special soundings ? 

I Input special surface obs ? 

! Number of stations only to be used 
! Station ID’S used 


1 Number of stations to be excluded 
! Station ID'S to be excluded 
» Prefix with 'r' for rawindsonde, 


’s' for surface 


! Minimum surface station separation in degrees. 

! Any surface obs within this distance ^ 

! of another obs will be thrown out unless it has 

! less missing data, in which case the other obs 
I will be thrown out. 

! If ISTAPLT = 1, soundings are plotted; 
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ISTAREP 

= 

0, 

i If ISTAREP = 1, soundings are listed; 





I no objective analysis is done. 





! If ISTAREP/ I STAPLT = 0, normal processing is done 

lUPPER 

= 

1. 

! 1-Do upper air analysis, O-Only do surface 

IGRIDFL 

- 

1, 

! Grid flag =0, if no grid point data. 

only obs 




! 1, if all grid point data 

and obs 




! 2, if partial grid point and obs 




! 3, if only grid data 


GRIDWT 

_ 

.01 

, .01, .01, 





! Relative weight for the gridded pressure data compared 




! to the observational data in the objective analysis 

GOBSEP 


5.0, 

1 Grid -observation separation (degrees) 


GOBRAD 

= 

5.0, 

! Grid- obs proximity radius (degrees) 


WVLNTH 


1600. , 

1200., 400., 





! Used in S. Barnes objective analysis. 





I Wavelength in km to be retained to 

the RESPON 




1 percent from the data to the upper 

air grids. 

SWVLNTH 

= 

1000. , 

400., 100., 





! Wavelength for surface objective analysis 

RES PON 


.90. 

.9, .9, 





! Percentage of amplitude to be retained. 

$END 





Graphical processing 


$ISAN_GRAPH 





1 

Main switches for plotting 


IPLTPRS 

— 

0, 

1 Pressure coordinate horizontal plots 


IPLTISN 

= 

1, 

1 Isentropic coordinate horizontal plots 


IPLTSIG 

= 

1, 

' Sigma- z coordinate horizontal plots 


IPLTSTA 

= 

0. 

! Isentropic coordinate "station" plots 



J 

Pressure 

plotting information 


ILFTII 

— 

0, 

1 Left boundary window 


IRGTII 

= 

18. 

i Right boundary window 


IBOTIJ 


• 3, 

! Bottom boundary window 


ITOPIJ 

= 

13, 

! Top boundary window 





! Window defaults to entire domain if one 

equals 0. 

NPLEV 


8. 

1 Number of pressure levels to plot 


IPLEV 

= 

1000,850, 

700,500,400,300,200,100, 





! Levels to be plotted 


NFLDUl 

- 

4, 

i Number of fields to be plotted 


IFLDUl 

= 

•U' , ’THETA' . 'GEO' RELHUM' , 1 Field names 


CONUl 

= 

0.,0.,0., 

0., ! Field contour increment 


IVELUl 


o 

o 

o 

1 Velocity vector flag 



1 

Isentropic plotting information 


ILFT3I 

— 

0, 

! Left boundary window 


IRGT3I 

= 

18, 

! Right boundary window 


IBOT3J 

= 

3, 

! Bottom boundary window 


ITOP3J 

= 

13, 

i Top boundary window 





! Window defaults to entire domain if one 

equals 0. 




! Upper air plots: 


IUP3BEG 

= 

340, 

! Starting isentropic level for plotting 


IUP3END 

= 

400, 

! Ending isentropic level 


IUP3INC 

= 

60, 

1 Level increment 


NFLDU3 

= 

4, 

! Number of fields to be plotted 
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IFLDU3 

C0NU3 

IVELU3 


•U' , ’ PRESS 

0 . , 0 . , 

1 , 0 , 


'MONTSTR* , 'RELHUM' , ! Field names 

I Field contour increment 
1 Velocity vector flag 


Surface plotting information 


Uses isentropic plotting window info 

4 1 Number of surface fields to plot 

'U' , 'TEMP' PRESS' , 'RELHUM' , i Fieldnames 

0.,0w0.,0.,0., ! Field contour increment 

l’o,o'o,'6,' I Velocity vector flag 


NFLDS3 = 
IFLDS3 = 
CONS3 
IVELS3 == 


! Sigma -z plotting information 

I - 

Uses isentropic plotting window info 


ISZBEG = 2, 
ISZEND = 26 
ISZINC = 6, 


! Starting sigma -z level for plotting 
! Ending sigma -z level 
! Level increment 


NFLDSZ 

IFLDSZ 

CONSZ 

IVELSZ 


4, 
'U' 
0 . , 0 . 
1 , 0 , 


I Number of fields to be plotted 
PRESS' , 'THETA' , 'RELHUM' , i Field names 

! Field contour increment 
! Velocity vector flag 


"Station" plotting information 


NPLTRAW=0, 


NSTIS3 

ISTIS3 


4, 

'U' , ' PRESS 


! number of station surface plots 
■ RELHUM TEMP ' , ! fieldnames 


Cross ' sec tion plotting information 


NCROSS3 


0, 

ICRTYP3 

= 

2,1, 

ICRA3 

= 

1,1, 

ICRB3 


35,43, 

ICRL3 

= 

22,25, 

NCRFLD3 

= 

3, 

ICRFLD3 

= 

'MIXRAT ' , ' RELHUM' 

THCON3 

= 

5. .5. , 5. , 

ACON3 

= 

0 . , 0 . , 0 . , 


! number of cross section slabs 
1 type of slab: 1=E-W, 2=N-S 
! left window 
! right window 
! cross section location 
! number of plots on each cross 
'THETAE', ! field names 

I contour interval of 
I contour interval of 


$END 

1 

! Field values for graphical stage 


Pressure Isentropic Station 


Sigma - z 


U 

U 

V 

V 

TEMP 

PRESS 

GEO 

GEO 

RELHUM 

RELHUM 

MIXRAT 

MIXRAT 

THETA 

THETA 

SPEED 

SPEED 

ENERGY 

ENERGY 

THETAE 

THETAE 

SPRESS 

SPRESS 


U U 

V V 

PRESS PRESS 

temp THETA 

RELHUM RELHUM 

MIXRAT 


section 

i sen tropes 
other field 
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APPENDIX B 
PROWESS RAMSIN file 
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&PARA 

IPARALLEL=1. 


NMACH=6 , HOSTS= ' par7 ' , • par6 ' , ' par5 ’ . ' par4 ' 
HPERF=1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 
HPERF=1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 
CPROG= ' /parO/ul/prowess/ system/bin/4an ' , 

CNAME=‘/parO/ul/prowess/system/fore/RAMSIN.op' 

CVTAB= ' /parO/ul/prowess/ system/ f ore/VTABLE ’ , 


par3 ' , 'par2 ’ , 'pari ’ , 


/ $END 


$MODEL_GRIDS 

1 Simulation title (64 chars) 

EXPNME = ’Version 3b sample' , 


lOTYPE 

0, 


j - i=diagnos t ic , 0-normal run 

RUNTY PE = 

* INITIAL' . 

I type of run: 

MAKESFC, INITIAL, HISTORY, or MAKEVFILE 

TIMEUNIT = 

'h' , 

t 

'h*,'m',*s' - Time units of 

TIMMAX, TIMSTR, VTIME 

TIMMAX= 

20.0, 

j 

Final time of simulation 

NGRIDS = 

4, 


! Number of grids to run 
1 

NNXP 

NNYP = 

NNZP 

34, 34, 38, 

30, 38, 38, 

30, 30, 30, 

42, 

58, 

30, 

I Number of x gridpoints 
! Number of y gridpoints 
1 Number of z gridpoints 

NNZG 

11,11,11,11, 


! Number of soil gridpoints 

NXTNEST = 

1, 1, 2, 3, 


! Grid number which is the next 
! coarser grid 


! Coarse grid specifications 


IHTRAN = 
DELTAX = 
DELTAY = 

DELTAZ = 

DZRAT 

DZMAX 


1 , 

72000. , 
72000. , 

75. , 
1.25, 
1000 . , 


! 0 -Cartesian, 1- Polar stereo 

! X and Y grid spacing 

! Z grid spacing (set to 0 . to use ZZ) 

! vertical grid stretch ratio 
! maximum delta Z for vertical stretch 


ZZ=0., 25., 75., 150., 250., 500. 
3250., 4000., 5000., 6000., 


750., 1000., 1500., 2000., 2500., 
! Vertical levels if DELTAZ = 


0 


DTLONG = 90., 

NRATIO = 4 , 


! Coarse grid long timestep 
1 Small timestep ratio 


IMONTHl 

IDATEl 

lYEARl 

STRTIM 


NSTRATX 

NSTRATY 

NNDTRAT 


01 , 

17, 


! Month 
Day 
Year 



12 

0, 

1 

GMT of model TIME = 0 . 






Nest ratios between this grid 






and the next coarser grid. 

_ 

1, 

4 , 

3, 4, 


x-direction 

_ 

1,* 

4, 

3, 4, 


y- direction 


1, 

2, 

2, 2, 


time 


NESTZ =0, 

NSTRATZ=4 ,4, 3, 3, 3, 3, 2, 2, 2, 2, 1,1, 1,1, 


contort coarser grids if negative 


POLELAT = 30.0, 

POLELON = -83.0, 

CENTLAT= 30.0,28.5,28.5,28.6, 
CENTLON= -83.0, -82.0, -81.0, -80.75, 


! Latitude of pole point 
! Longitude of pole point 


! Grid point on the next coarser 
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! nest where the lower southwest 
I corner of this nest will start. 

1 IF NINEST or NJNEST = 0, use CENTLAT/LON 


NINEST 

= 

1, 

0, 

0, 0, 


I i- point 

NJNEST 

= 

1, 

0, 

0. 0, 


! j -point 

NKNEST 

= 

1, 

1, 

1, 1. 


I k-point 

NNSTTOP 


1, 

1 

1, 

1, 

! Flag (0-no or 1-yes) if this 

NNSTBOT 


1, 

1, 

1. 

1, 

1 nest goes the top or bottom of 
I coarsest nest. 

GRIDU 

= 

0. , 

0. , 

, 0. , 

0. , 

1 u- component for moving grids 

GRIDV 


0. , 

0. , 

, 0. , 

0. , 

J v- component for moving grids 


(still not working*) 


$END 


$ MODEL_F I LE_I NFO 


! History file input 

TIMSTR= 1500., ! time of history start (see TIMEUNIT) 

HFILIN = ’ cg5 .hlSOOs ' , ! input history file name 


! History/analysis file output 


IOUTPUT= 1, ! 0-no 

HFILOUT='hist/f .h‘ , 

AFILOUT= ' f . a ' , 

HFUNITS= ' (3 ' , AFUNITS= ' d ’ , 

FRQHIS =72000., FRQANL =3600., 

! Variable initializ 


files, 1-save in ASCII, 2 -save in binary 
! history file prefix 
! analysis file prefix 
! history/anal file units (M,m, H, h, S, s) 

! history/anal file frequency 

tion input 


INITIAL = 2, 


I Initial fields - l=horiz , homogeneous , 2=variable 


VARFIL(l) =' isan/pvl7 - jan-97 -00 ' , 
VARFIL ( 2 ) = ' isan/pvl7 - j an - 97 - 06 • , 
VARFIL ( 3 ) = ' i san/pvl 7 - j an - 9 7 - 12 ' , 
VARFIL (4 )=’isan/pvl7- jan-97 -18 ' , 
VARFIL (5) = ' i san/pvl 7 - jan-97 - 24 ' , 


VTIME(l) =00. , 
VTIME{2) =06. , 
VTIME(3)=12. , 
VTIME(4) =18. , 
VTIME(5)=24. , 


IVWIND = 0, 


» initial winds ( only 0 works) 


NUDLAT 

TNUDLAT 

TNUDCENT 

TNUDTOP 

ZNUDTOP 


5, 

900. , 

0 . , 

0 . , 

15000. , 


I number of points in the lateral bnd region 
I nudging time scale (s) at lateral boundary 
nudging time scale (s) in center of domain 
nudging time scale (s) at top of domain 

! nudging at top of domain above this height (m) 


! Printed output controls 


FRQPRT = 21600., 
FRQIPR = 99930., 
FRQIST = 99900., 


Printout frequency 
Integral print frequency 
Integral store frequency 


ISTPFL 

INITFLD 

INPRTFL 


1, ! Timestep message frequency flag 

0, ! Initial field print flag 0=no print, l=print 

1, ! Namelist print flag 0=no print, l=print 


! Input topography variables 


SFCFILES = ' /parO/ul/prowess/system/f ore/sf c2 ' , 
surface characteristic 


! files. 


File path and prefix for 


ITOPTFLG 

IPCTLFLG 

ISSTFLG 

IVEGTFLG 


1 . 1 , 1 , 1 , 
1 . 1 , 1 , 1 , 
1 . 0 , 0 , 0 , 
1 , 1 , 1 . 1 , 


2 - Fill data in ’*ruser“ 

1 - Interp data from latlon dataset 
0 - Interpolate from coarser grid 


The following only apply for IxxxxFLG=l 
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ITOPTFN 


I PCTLFN 


' /ul/prow6SS / sys t6ir\/ data./ rams/ topolOm/H 
' /ul/prowess/sys tern/ data/ rams/ topo30s/U 
' /ul/prowess/ system/ data/rams/ topo30s/U 

' /ul/prowess/ sys tem/data/rams/topo30s/U 

» /ul/prowess/ system/ data/ rams/pctllOm/L 

• /ul/prowess/ system/data/rams /pctllOm/L 

* /ul /prowess/ sys tem/data/ rams/pctllOm/L 

1 /ul/prowess/ sys tern/ da ta/rams/pctllOm/L 

ISSTFN = ' /ul /prowess/ sys tern/ da ta/rams/sst/ SJUL , 

IVEGTFN = ’ /ul/prowess/system/data/rams/ndvi/V , 

> /ul /prowess/ sys tem/data/rams/ndvi/V , 

* /ul/prowess/sys tem/data/rams/ndvi/V , 

• /ul/prowess/sys tem/data/rams/ndvi/V , 


Input topography file name 
Input topography file name 
Input topography file name 


SILAVWT 

TOPTWVL 

PCTLWVL 

SSTWVL 


0.0, 

0. , 

, 0. ,0 

4.0, 

4. , 

, 4. 

,4 

2.0, 

2. , 

, 2. 

,2 

2.0, 

2. , 

r 2. 

,2 


Weighting of topo silhouette averaging 
Topo wavelength cutoff in filter 
Land pet wavelength cutoff in filter 
Land pet wavelength cutoff in filter 


MKCOLTAB =0, I make microphysics collection table: 0 

COLTABFN = ' /par0/ul/prowess/system/rams3b/data/micro/ct2 . 0 , 

filename to read or write ok 0» 

EVPTABFN = Vpar0/ul/prowess/system/rams3b/data/micro/et2.u , 

filename to read 


no, 1 = yes 
! collection table 

! evaporation table 


$END 


$MODEL_OPTIONS 

NADDSC 


0 , 0 , 0 , 0 , 


Number of additional scalar species {< 5) 


NTOPSMTH = 
I Z FLAT = 


0 , 

0 , 


Number of passes for topography smoother . . 

width of flat margin around domain {in grid points) 


Numerical schemes 


ITMDIFF 

= 

3, 

NONHYD 

= 

1, 

SSPCT 


0. , 

IMPL 

= 

1, 

ICNTEQ 

= 

2, 

WTKD 

= 

0. , 

ICORFLG 


1, 

IBCTOP = 


0, 

IBND 

_ 

1, 

JBND 


1, 

CPHAS 

= 

20. , 

LSFLG 


0, 

NFPT 


0, 

DISTIM 

— 

60. , 

TIMSCL 


0. , 

KSPIN 

= 

33, 

KMSPIN 

= 

10, 

IPRSPLT 


0, 

lADVL 

= 

2, 

lADVF 


2, 

IPGRAD 

= 

1, 

FILT4 


000. 

FXLONG 

= 

0.0000 


l=forward, 2=leapfrog, 3-hybrid 
nonhydrostatic=l, hydrostatic=0 

Sound speed fraction for the nonhydros ta tic model 
Implicit flag for acoustic model - 0=off, 

Hyd - continuity equation - 1 -incomp, 2 -anelastic 
Klemp/Durran - current times tep weighting 

Coriolis flag/2D v-component - 0=off, l=on 

top boundary condition ^ 

0- wall on top(nonhyd) 2 -prognostic sfc prs(hyd) 

1 - Klemp-Durran (nonhyd/hyd) 3 -material surface (hy 

Lateral boundary condition flags 
1 -Klemp/Wilhelmson, 2 -Klemp/Lilly, 3-Orlanski 
4 - cyclic 

Phase speed if IBND or JBND =1 . 

Large-scale gradient flag for variables other than 
normal velocity: 

0 = zero gradient inflow and outflow 

1 = zero gradient inflow, radiative b.c. outflow 

2 = constant inflow, radiative b.c. outflow 

3 = constant inflow and outflow 

Rayleigh friction - number of points from the top 
- dissipation time scale 


Initial wind spin- up 


time scale 

below this level {velocity) 
- below this level (momentum) 

precipitation time-split scheme - 0=°“' 

Order of advection - Leapfrog - 2 or 4 
Order of advection - Forward - 2 or 6 
Pressure gradient scheme for topo.: 

Fourth order filter coefficient (0 - off) 

{ > 0 times tep at which 2 delta waves are 
totally removed) 

Long filter coefficient 0=off 
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FYLONG = 0.0000, I Long filter coefficient 0=off 

! Radiation parcuneters 

ISWRTYP = 1, ! Shortwave radiation type 0-none, 2 -Mahrer/Pielke, 1-Chen 

ILWRTYP = 1, ! Longwave radiation type 0-none, 2 -Mahrer/Pielke, 1-Chen 

RADFRQ = 1200., ! Frequency of radiation tendency update in seconds 

LONRAD = 1, • Longitudinal variation of shortwave (0-no, 1-yes) 

! Cumulus parameterization parameters 

NNQPARM =1, 1, 1, 0, 5 convective parameterization flag 

! 0-off, 1-on 

CONFRQ = 1200., I Frequency of conv param. updates in seconds 

WCLDBS = .001, ! vertical motion needed at cloud base for convection 

! Surface layer and soil parameterization 

ISFCL = 1, ! surface layer/soil/veg model 

! 0- specified surface layer gradients 
I 1 - soil/vegetation model 

NVGCON = 1, ! Vegetation type 


1 -- 

Crop/mixed farming 

2 -- 

Short grass 

3 - - 

Evergreen needleleaf tree 

4 - - 

Deciduous needleleaf tree 

5 - - 

Deciduous broadleaf tree 

6 - - 

Evergreen broadleaf tree 

7 - - 

Tall grass 

8 - - 

Desert 

9 - - 

Tundra 

10 - - 

Irrigated crop 

11 -- 

Semi - desert 

12 -- 

Ice cap/glacier 

13 -- 

Bog or marsh 

14 -- 

Inland water 

15 -- 

Ocean 

16 - ■ 

Evergreen shrub 

17 - - 

Deciduous shrub 

18 -- 

Mixed woodland 


TSEASN = 298., ! Average seasonal temp 

TVGOFF = 0., ! Initial veg temp offset 

VWTRCON = -000, ! Initial veg water storage (not working) 

UBMIN = 0.25, i Minimum U value to use in computing 

PCTLCON =1.0, ! constant land percentage if for all domain 

NSLCON =6, ! constant soil type if for all domain 

! l=sand 2=loamy sand 3=sandy loam 

! 4=silt loam 5=loam 6=sandy clay loam 

I 7=silty clay loam 8=clay loam 9=sandy clay 

I 10=silty clay ll=clay I2=peat 

ZROUGH =0.05, ! constant roughness if for all domain 

ALBEDO = 0.2, ! constant albedo when not running soil model 

SEATMP = 280., 1 constant water surface temperature 

DTHCON = -0., ! constant surface layer temp gradient for no soil 

DRTCON = .000, ! constant surface layer moist gradient for no soil 

SOILDZ =0., I soil model grid spacing 

SLZ = - .50, - .40, - .30, - .25, - .20, - .16, - .12, - .09, - .06, - .03, 0.,! soil grid levels 

SLMSTR = 0.35, 0.35, 0.35, 0.35, 0.35, ! initial soil moisture 

0.35, 0.35, 0.35, 0.35, 0.35, 0.35, 

STGOFF= 5,, 5,, 5., 5., 3.5, 2., .5, - 1 . , - 1 . 5 , - 1 . 8 , -2., 

! Initial soil temperature offset from lowest 
! atmospheric level 


! Eddy diffusion coefficient parameters 

IDIFFK = 1,1, 1,1, IK flag: 

! 1 - Horizontal deformation/ Vertical Mellor - Yamada 

1 2 - Anisotropic deformation (horiz and vert differ) 

! 3 - Isotropic deformation (horiz and vert same) 

! 4 - Deardorff TKE (horiz and vert same) 

CSX = .20, ,20, .20, .20, .20, ! Adjustable parameter, deformation horiz. K's coefficient 

CSZ = . 20, . 20, . 20, . 20, . 20 , ! Adjustable parameter, deformation vert- K’s coefficient 

XKHKM = 3. ,3. ,3. ,3. ,3., ! Ratio of horizontal K_h to K_jn for deformation 

ZKHKM = 3. ,3. ,3, ,3. ,3., ! Ratio of vertical K_h to K_m for deformation 

AKMIN= 1.0, .75,1.0,2.0,1., I Ratio of minimum horizontal eddy viscosity coefficient 

! to typical value from deformation K 

! Microphysics 

NLEVEL = 3, 3, 3,3, 1 moisture complexity level 


NLEVEL = 

1, 

1, 

1,1, 

I moisture complexity level 

INUCPRG = 

0, 

0, 

o 

o 

! Prognose ice nuclei (0-no, 1-yes) 

ICLOUD = 

0, 

0, 

o 

o 

! Microphysics flags 
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I RAIN 

= 

0, 

0, 0,0, 


IPRIS 

= 

0, 

0, 0,0, 


I SNOW 

= 

0, 

0, 0,0, 


lAGGR 

= 

0, 

0, 0,0, 


IGRAUP 

= 

0, 

0, 0,0, 


I HAIL 

= 

0, 

0, 0,0, 


I CLOUD 

- 

4, 

4, 4,4, 

1 

IRAIN 

= 

2, 

2, 2,2, 

! - 

IPRIS 

- 

5, 

5, 5,5, 

1 

I SNOW 

= 

0, 

0, 0,0, 

1 

lAGGR 

= 

2, 

2, 2,2, 

t 

IGRAUP 

= 

0, 

0, 0,0, 

j 

I HAIL 

= 

2, 

2, 2,2, 

j 

1 

CP ARM 

= 

,3e9, 

.3e9, .3e9, 

.3e9, 

RPARM 

= , 

.le-2, 

,le-2, .le-2, 

. le- 2 , 

PPARM 

= 

0. , 

o 

o 

0. , 

SPARM 

= 

.le-2. 

.le-2, .le-2. 

.le-2. 

APARM 


. le-2. 

.le-2, -le-2. 

. le- 2 , 

GPARM 

- 

. le- 2 , 

.le-2, .le-2, 

.le-2. 

HPARM 

= 

•3e-2, 

.3e-2, .3e-2, 

. 3e - 2 , 

AMIO 

= 

1 .e- 

12, 

! mir 

gnus for; 
GNU = 

cloud 
2. , 

rain pris snow aggr 
2. , 2. , 2. , 2. , 


0 - no species 

1 - diagnostic concen. 

2 - specified mean diameter 

3 - specified y- intercept 

4 * specified concentration 
Microphysics flags 


0 - no species 

1 - diagnostic concen. 

2 ' specified mean diameter 

3 - specified y- intercept 

4 . specified concentration 

5 - prognostic concentration 

Microphysics parconeters 


Characteristic diameter, 

number concentration, or 
y- intercept 


2 . , 2 . 


gamma shape parms 


$END 

$MODEL_SOUND 

! Sounding specification 


IPSFLG=1, 
ITSFLG=2 , 
IRTSFLG=3, 

lUSFLG^O, 
lUSRC = 0, 


Flags for how sounding is specified 

! specifies what is in PS array 
! 0 " pressure (mb) , 1 ■ heights (m) 


PS(l)=sfc press (mb) 


specifies what is in TS array 
0 - temp(C), 1 - temp(K), 2 - pot. temp(K) 


' specifies what is in RTS array , \ 
! 0 - dew pnt.(C), 1 - dew pnt.(K), 2 - mix rat(g/kg) 
I 3 - relative humidity on %, 4 


dew pnt depression (K) 


* specifies what is in US and VS arrays 

1 0 - u,v component (m/s) , 1 - umoms -direction, vmoms-speed 

! source of wind profile: 

! 0 - umoms, vmoms, valid at sounding levels (PS) 

» -1 - usndg, vsndg, valid at model levels (Z) 

HS = 0. , 

PS = 1000., 1000. ,2000. ,3000. ,3500. ,4500. ,7000., 10000. ,15000., 20000. ,25000. 


TS = 

300. 

,302. ,304. ,306. ,308 

. ,310. 

,314. , 320. ,370. ,420. ,470., 

RTS = 

60. 

,60. , 60 

. , 60 . , 10 . , 10 

. , 10. , 

10. , 10 . , 10 . , 10. , 

US = 

3. ,3 

. ,3. ,3. 

,3. ,3. ,3. ,3. 

,3. ,3. 

,3. , 

VS = 

0. ,0 

. ,0. ,0. 

o 

o 

o 

o 

,0. ,0. 

. ,0. , 

USNDG 

VSNDG 

II II 
o o 

, 1 Wind components if 

t 

IUSRC= “1 

KMEANl =0, ! 

KMEAN2 =0, ! 

lower model 
upper model 

level 

level 

for calculation of (umean 
for calculation of (umean 

UMEAN 

VMEAN 

= 0. 
= 0. 

o o 

u - component 
V - component 

for Galilean transformation 
for Galilean transformation 


$END 


$MODEL PRINT fields to be printed during the simulation 

nplt = 6 ! Number of fields to be printed at each time 

! for various cross-sections (limit of 50) 
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IPLFLD 


= ' UP ' , • VP ' , • WP • . ■ THETA • , ’ TKE ’ , ’ RT ' , 

! Field names • see table below 
!! PLFMTd) = '10PF7.2', ! Format specification if default is unacceptable 

IXSCTN = 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 

! IXSCTN = 2, 2, 2, 2, 2, 2, 2.2,2, 

I Cross-section type (1=XZ, 2=YZ, 3=XY) 

ISBVAL = 2, 2, 2, 2, 2, 2, 2, 2. 2, 2, 2, 

! Grid -point slab value for third direction 

! The following variables can also be set in the $PRNT namelist: lAA, 

! lAB, JOA, JOB, NAAVG, NOAVG, PLTIT, PLCONLO, PLCONHI, and PLCONIN. 


$END 


C 

•UP' 

- UP(M/S) 

'RC 

- RC(G/KG) 

' PCPT' 

- TOTPRE 

c 

1 VP' 

- VP(M/S) 

'RR' 

■ RR(G/KG) 

•TKE' 

- TKE 

C 

'WP' 

- WP(CM/S) 

'Rp. 

■ RP(G/KG) 

•HSCL' 

- HL(M) 

C 

' PP' 

PRS (MB) 

'RA' 

RA(G/KG) 

'VSCL' ■ 

- VL(M) 

c 

•THP’ 

THP(K) 




TG (K) 

C 

’ THETA ' 

THETA (K) 

'RL' 

RL(G/KG) 

'TG' 

c 

•THVP' 

THV (K) 

•RI ■ 

RI (G/KG) 

■ SLM’ 

SLM (pct; 

c 

•TV 

TV(K) 

’ RCOND • - 

RD(G/KG) 

' CONPR ' 

CON RATE 

c 

•RT’ 

RT(G/KG) 

•CP' 

NPRIS 

■ CONP • 

CON PCP 

c 

'RV 

RV(G/KG) 

,RTp. 

RT’ (G/KG) 

' CONH • 

CON HEAT 

c 





'CONM' 

CON MO IS 

c 

'THIL' 

- Theta- il (K) 

'TEMP 

' - temperature (K) 



Q ‘TVP* - Tv‘ (K) 'THV' - Theta -V (K) 

C ' RELHUM’ - relative humidity {%} 'SPEED' - wind speed (m/s) 

C 'FTHRD'- radiative flux convergence (] 

C 'MICRO'- GAS PRC 

C 'ZO' - ZO (M) 'ZI' - ZI (M) 'ZMAT' - ZMAT (M) 

C ' USTARL ' -USTARL (M/S) ' USTARW* - USTARW (M/S) ' TSTARL ' -TSTARL (K) 

C ' TSTARW -TSTARW(K) ' RSTARL ' - RSTARL (G/G) ' RSTARW - RSTARW (G/G) 

C 'UW' - UW (M*M/S*S) 'VW' - VW (M*M/S^S) 

C 'WFZ' - WFZ (M*M/S*S) ' TFZ ’ - TFZ (K*M/S) 

C 'QFZ' - QFZ (G*M/G*S) ’ RLONG ' - RLONG 

C ' R SHORT ' -RSHORT 


$ISAN_CONTROL 

MSTAGE = 1, 1, 1, ! Main switches for 

! pressure, isentropic. "varfile" processing 

NATIMES =1, 1 Number of times on which to perform analysis 

lAHOUR =24, ! Hours to analyze 

lADATE =17, ! Dates 

lAMONTH =01, i Months 
lAYEAR =97, ! Years 

GUESS1ST=' PRESS' , 

! NCAR archived data file names 

lAPR = ' . . /data/ipl7- jan-97 -24 ' , 

lARAWI = ' ' , ! Archived rawindsonde file name 

lASRFCE = ’ ' , 1 Archived surface obs file name 

J File names and dispose flags 
IFNPRS = 'isan/pp', ! Pressure file name prefix 

IOFLGP= 1, ! Dispose flag: 0 = no write, 1 = write 

IFNISN = 'isan/pi', ! Isentropic file name prefix 

IOFLGI= 1, ! Dispose flag: 0 = no write, 1 = write 

IFNSIG = "isan/ps’, I Sigma-z file name prefix 

I0FLGS= 1, I Dispose flag: 0 = no write, 1 = write 

IFNVAR = 'isan/pv, ! "Variable initialization" file name prefix 
I0FLGV= 1, ! Dispose flag: 0 = no write, 1 = write 


$END 


80 


$ISAN_PRESSURE 


NPRX=93, 

NPRY=65, 

NPRZ=19, 


Pressure grid information: 

] number of grid points in x (Ion) direction 
! number of grid points in y (lat) direction 
! number of pressure levels 


WPLON = -140.0, 

SPLAT = 20.0, 


I West longitude bound of pressure data access 
! South latitude bound , j *. 

! Latitude and longitude bounds are north and east 
! if positive, south and west if negative. 


SPCNPRX = 2.5, 
SPCNPRY = 1.25, 


! East -west grid spacing of pressure data 
! North -south spacing 


LEV PR 


! Pressure levels (mb) in input dataset to access 
1000,950 ,900,850,800,750,700,650,600, 

550,500,450,400, 350, 300,250,200, 150, 100, 


$END 


Isentropic and sigma- z processing 


$ I SAN_I SENTROPIC 


Specify isentropic levels 


NISN 

LEVTH 


32 ! number of isentropic levels 

30,282,284,286,288,290,292,294,296,298,300,303,306,309,312, 
15,318,321,324,327,330,335,340,345,350,355,360,380,400,420, 
n A c 


Analyzed grid information: 


NIGRIDS = 3, 

TOPSIGZ = 20000 

HYBBOT = 4000., 
HYBTOP = 6000., 

SFCINF = 1000., 

SIGZWT = 1 . , 


I number of RAMS grids to analyze 
1 sigma- z coordinates to about this height 

! Bottom (m) of blended sigma - z/isentropic layer in varfiles 
! Top (m) of blended sigma - z/isentropic layer 

I Vertical influence of surface observation analysis 


Weight for sigma -z data in varfile; HvnnnT 

n = no siaz data, l.=full weight from surface to HYBBOT 


NFEEDVAR = 1, 


MAXSTA=150, 
MAXSFC=1000 , 

lARCSND = 1, 

lARCSFC = 1 , 

ISPCSND = 0, 

ISPCSFC = 0, 

NONLYS = 0 , 

IDONLYS =' 76458 ' , 

NOTSTA = 0 , 

NOTID ='r76458', 


STASEP = .001, 


I 1 = feed back nested grid varfile info, 0 - don’t 


Observation number limits: 


! maximum number of rawindsondes (archived special) 
! maximum number of surface observations 

! Input archived soundings ? (0-no, 1-yes) 

! Input archived surface obs ? 

! Input special soundings ? 

I Input special surface obs ? 

I Number of stations only to be used 
I Station ID's used 

! Number of stations to be excluded 
i Station ID's to be excluded 

« Prefix with *r’ for rawindsonde, 's' for surface 


! Minimum surface station separation in degrees. 
Any surface obs within this distance 
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of another obs will be thrown out unless it has 
less missing data, in which case the other obs 
will be thrown out. 


ISTAPLT 


0, 


I If ISTAPLT = 1, soundings are plotted; 

I STARE P 


0, 


! If ISTAREP = 1, soundings are listed; 

! no objective analysis is done. 

! If ISTAREP/ ISTAPLT = 0, normal processing is done 

lUPPER 

- 

1, 


! 1-Do upper air analysis, O-Only do surface 

IGRIDFL 

_ 

1, 


I Grid flag =0, if no grid point data, only obs 

IGRIDFL 


3. 

! Grid flag =0, if no grid point data, only obs 
t 1, if all grid point data and obs 

! 2, if partial grid point and obs 

j 3, if only grid data 

GRIDWT 


.01, .01 

, .01, 

! Relative weight for the gridded pressure data compared 
' to the observational data in the objective analysis 

GOBSEP 


5.0, 


1 Grid- observation separation (degrees) 

GOB RAD 

= 

5.0, 


! Grid- obs proximity radius (degrees) 

WVLNTH 

“ 

1200. , 

900. , 

600., 

! Used in S. Barnes objective analysis. 

! Wavelength in km to be retained to the RESPON 

• percent from the data to the upper air grids. 

SWVLNTH 


750. , 

300. , 

150. , 

! Wavelength for surface objective analysis 

RES PON 


.90, 

.9, 

.9, 

! Percentage of amplitude to be retained. 


$END 


Graphical processing 


$ISAN„GRAPH 

! Main switches for plotting 


IPLTPRS 

= 

0, 

Pressure coordinate horizontal plots 


IPLTISN 

= 

1, 

Isentropic coordinate horizontal plots 


IPLTSIG 

= 

1, 

Sigma -z coordinate horizontal plots 


IPLTSTA 

- 

0, 

Isentropic coordinate ’’station" plots 



1 

Pressure plotting information 


ILFTII 


0. ! 

Left boundary window 


IRGTII 

= 

18, ! 

Right boundary window 


IBOTIJ 


3, ! 

Bottom boundary window 


ITOPIJ 

= 

13, ! 

Top boundary window 

equals 0 



1 

Window defaults to entire domain if one 

NPLEV 


2, 

i Number of pressure levels to plot 


IPLEV 

= 

1000,500, 

! Levels to be plotted 


NFLDUl 

= 

4, 

! Number of fields to be plotted 


IFLDUl 

= 

• U • , ' THETA ’ , ' GEO ' , ' RELHUM • , ! Field names 


CONUl 

= 

0. ,0. ,0. ,0 

, [ Field contour increment 


IVELUl 


2, 0,0,0, 

1 Velocity vector flag 



j 

I sen tropic 

plotting information 


ILFT3I 

1 

0, ! 

Left boundary window 


IRGT3I 

- 

18, • 

Right boundary window 


IBOT3J 


3, I 

Bottom boundary window 


ITOP3J 

= 

13, 1 

Top boundary window 

equals 0 


1 

Window defaults to entire domain if one 



r 

Upper air plots: 
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IUP3BEG = 340, 1 Starting isentropic level for plotting 

IUP3END = 4000, 1 Ending isentropic level 

IUP3INC =60, 1 Level increment 

MFT DU3 = 4 ‘ Number of fields to be plotted 

IFLDU3 = -U* , ’PRESS' , ’MONTSTR' , -RELHUM’ , I Fieldnames 

CONU3 =0,0., ^ Field contour increment 

IVELU3 = l!o,*’ I Velocity vector flag 


Surface plotting information 


1 Uses isentropic plotting window info 

NFTns3 = 4 I Number of surface fields to plot 

IFLDS3 = ’U' , 'TEMP- , 'PRESS’ , 'RELHUM' , 1 Fieldnames 

CONS3 = 0 ., 0 ., 0 ., 0 ., 0 ., ! Field contour increment 

IVELS3 = 1 , 0 , 0 , 0 , 0 , ! Velocity vector flag 


! Sigma 'Z plotting information 

I 

! Uses isentropic plotting window info 


ISZBEG = 2, 
ISZEND = 26, 
ISZINC = 6, 

NFLDSZ = 4, 

IFLDSZ = ’U’ 

CONSZ = 0 . , 0 . 
IVELSZ = 1,0, 


! Starting sigma- z level for plotting 
I Ending sigma- z level 
! Level increment 

! Number of fields to be plotted 
'PRESS’ , 'THETA’ , ’RELHUM’ , i Field names 

! Field contour increment 
I Velocity vector flag 


NPLTRAW 


''Station" plotting information 


1 Approximate number of raw rawinsonde 
1 frame. 0 turns off plotting. 


plots per 


NSTIS3 

ISTIS3 


3 ^ ! number of station surface plots 

' PRESS' , 'TEMP' , ’RELHUM’ , ! fieldnames 


Cross-section plotting information 


NCROSS3 

= 

0, 

ICRTYP3 

= 

2,1, 

ICRA3 


1,1. 

ICRB3 

= 

35,43, 

ICRL3 


22,25, 

NCRFLD3 

= 

3, 

ICRFLD3 

= 

'MIXRAT ' , ' RELHUM 

THCON3 

= 

5. ,5. ,5. , 

ACON3 

= 

0. ,0. ,0. , 


! number of cross section slabs 
1 type of slab: 1=E-W, 2=N-S 
! left window 
! right window 
! cross section location 
I number of plots on each cross section 
'THETAE', ! field neimes 

! contour interval of isentropes 
I contour interval of other field 


$END 

I - - 

1 Field values for graphical stage 


Pressure 


Isentropic Station 


Sigma - z 


U 

U 

V 

V 

TEMP 

PRESS 

GEO 

GEO 

RELHUM 

RELHUM 

MIXRAT 

MIXRAT 

THETA 

THETA 

SPEED 

SPEED 

ENERGY 

ENERGY 

THETAE 

THETAE 

S PRESS 

SPRESS 


U U 

V V 

PRESS PRESS 

temp theta 

RELHUM RELHUM 

MIXRAT 
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APPENDIX C 
ERDAS HYPACT_IN file 
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&hypact_in 

hypref='/ul/met/9701171200.alt3/RAMS/f .a 
freqpv^3rt= 600.0000000 , 

hpartfile=' ./test/H. 
dtpart= 120.00000000 , 

inaxpart= 100000, 

npartbl= 10000, 

freqavg= 600.0000000 , 

npsources= 21, 

emission= 305.7910156 , 1242.040039 


emission=” 305.7910156 , 1242.040039 , 279.9200134 297.2650146 

641 5090332 , 695.7250366 , 629.8619995 , 655.7630615 , 675.3580322 

^182 3^0088 : ^350. 410156 , 672.8680420 , 637.9080200 , "98 '270264 , 

593 5370483 572.7090454 , 504.3070374 , 429.7100220 , 1379.640015 , 

1599 430054 3960 960205 , 0 . OOOOOOOOOOE+00 , 0 . OOOOOOOOOOE+00 , 0 . OOOOOOOOOOE+00 , 

o.ooooooooooE+00, o.ooooooooooE+oo, o.ooooooooooE+oo, o.ooooooooooE+oo, o-ooooooooooE+oo, 

0 . OOOOOOOOOOE+OO , O.OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO, O.OOOOOOOOOOE+00, 

0 OOOOOOOOOOE+OO 0 OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE 00, 
o;oZoZool+oo: o:ooooooooooE+oo: O.OOOOOOOOOOE+OO. O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 


695.7250366 

3350.410156 

572.7090454 

3960.960205 


O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

O.OOOOOOOOOOE+OO, 

s zpwr = 1.000000000 , 

szmin= 1.000000000 , 

szmax= 1.000000000 , 


O.OOOOOOOOOOE+OO, 


irplstrt= 163200, 163200, 163^UU, idjzuu, 

163200, 163200, 163200, 163200, 163200, 163200, 163200, 163200, 

163200, 163200, 163200, 163200, 163200, 163200, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 

0 0 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, ,,11 

ireldur= 1, 1. 1- i i 

1 1 1, 1. 1- 

r 1 i: 1. 1. 1' 

o! 0 , 0 , 0 , 0 . 0 , 0 , 0 , 

0 0 0, 0, 0, 0, 0, 0, 

0 . o! o! 0 , 0 , 0 , 0 , 0 , 

0, 0, 0, 

isimend= 203000, 

istrtdays= 0, 0, 0, 0, 0, 0, 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, A n I 

idurdays= 0, 0, 0, 'n 

0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0. 

0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 

ienddays= 0 , 

tricoord= O.OOOOOOOOOOE+OO, on c/toin-m 

srcx= -80.55101776 , -80.55413818 , -80.55320740 , -80.55086517 , -80.54810333 

, -80.54810333 , -80.54810333 , -80.54810333 , -80.54810333 -80.54810333 

80 54810333 , -80.54810333 , -80.54810333 , -80.54810333 , -80.54810333 , - 

80;54810333 , -80.54810333 , -80.54810333 , -80.54810333 , -80.54810333 - 

80 54810333 0 OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

0 OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

O.OOOOOOOOOOE+OO. O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O-OOOOOOOOOOE+00, 

0 OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

srcy= 28.42615509 , 28.42622375 , 28.42408943 , 28.42531013 , 28 42695618 

28.42696618 , 28.42695618 , 28.42695618 , 28.42695618 , 28 42695618 

28.42695618 . 28.42695618 . 28.42695618 , 28.42695618 , 28.42695618 

28 42695618 28.42695618 . 28.42695618 , 28.42695618 , 28.42695618 

28 ! 42695618 O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

0 OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, 

0 OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, O.OOOOOOOOOOE+OO, ^ ' 

srcz= 349 8999939 , 449.8999939 , 508.2500000 , 524.9499512 , 550.0000000 

583.4500122 , 614.4500122 , 643.0000000 , 671.5499878 , 742.9499512 


163200, 

163200, 

163200, 

0 , 


163200, 

163200, 

163200, 

0 , 


163200, 

163200, 

163200, 

0 , 


163200, 

163200, 

163200, 

0 , 


163200, 

163200, 

0 , 


1 , 

1 , 

0 , 

0 , 

0 , 

0 , 

203000, 
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856 3499756 . 925.0999756 , 950.0999756 , 975.0999756 , 1000.949951 

1028.599976 , 1057.199951 , 1085.750000 , 1200.000000 , 1450.099976 

2050.100098 , 0 . OOOOOOOOOOE+00 , 0 . OOOOOOOOOOE+00 » 0 . OOOOOOOOOOE+00 , 0 . OOOOOOOOOOE+00 , 

0 . OOOOOOOOOOE+00, 0. OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+00 , 0 . OOOOOOOOOOE+00 , 

0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO, 
0 . OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO, 
xsize= 470.3999939 , 705.7999878 , 789.0000000 , 808.2000122 ^833.79998' 


808.2000122 


833.7999878 


862 0000000 , 883.0000000 , 898.0000000 , 9U9.^UUUi^^ , yzi.zuuuxzz 

894.5999756 , 848.7999878 , 826.0000000 , 799,4000244 , 767.4000244 

727 5999756 , 679.4000244 , 622.2000122 , 399.7999878 , 399.7999878 

399!7999878 , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 , OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO, 0 .OOOOOOOOOOE+OO, 0 .OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO, 

0 . OOOOOOOOOOE+OO, 0 .OOOOOOOOOOE+OO, 0 .OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 , OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO, 

0 . OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0, OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 , OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 
ysize= 470.3999939 , 705.7999878 , 789.0000000 , 808.2000122 , 833.7999878 

RfiO nnnnnnn 883.0000000 , 898.0000000 , 909.2000122 , 921.2000122 , 


, 898.0000000 

826.0000000 
622.2000122 


, 909.2000122 

799,4000244 

399.7999878 


, 921.2000122 

767.4000244 

399.7999878 


0 . OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OC 
ysize= 470.3999939 , 705.7999878 , 789.000C 

, 862.0000000 , 883.0000000 , 898.0000000 

894.5999756 , 848.7999878 , 826.0000000 

727.5999756 , 679.4000244 , 622.2000122 

399.7999878 , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 


799.4000244 , 767.4000244 

399.7999878 , 399.7999878 

0. OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO, 


0 . OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 .OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0 .OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0. OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 , OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 
zsize= 100.0000000 , 100.0000000 , 16.69998169 , 16.70001221 , 33.40002441 


, 33.50000000 

112.5000000 

28.60003662 

899.8000488 

0. OOOOOOOOOOE+OO 


, 28.50000000 

25.00000000 

28.59997559 


, 28.59997559 

25.00000000 

28.50000000 


, 28,50000000 

25.00000000 

200.0000000 


, 114.2999878 

26.70001221 

300.1999512 


oyy.ouuu4oo , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 
0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 , OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 
0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 
0 . OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0 .OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO, 
0 . OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 .OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO, 


U . UUUUUUUUW , V , WUUUWU WV WWiJ ’ wvf , V , ^ ^ ^ ^ ^ ^ ^ ^ ^ . 

0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 0 , OOOOOOOOOOE+OO , 

0 . OOOOOOOOOOE+OO, 0. OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO , 0 . OOOOOOOOOOE+OO , 
__ _ ..... .. ^nnnnnr\nnr\t-%i/\f\ n r\r\Annnr\n A ncx nr\ 


0. OOOOOOOOOOE+OO, 0 . OOOOOOOOOOE+OO, o.c 
ihfall= 0, 

ihturb= 1, 

sourcetype ( 1) = * lag/ 1/ellipse/ total ’ , 
sourcetype (2 ) = ' lag/l/ellipse/ total ' , 
sourcetype (3 ) = ’ lag/l/ellipse/ total ' , 
sourcetype (4 ) = ’ lag/l/ellipse/ total ’ , 
sourcetype (5) = ’ lag/l/ellipse/ total ' , 
sourcetype ( 6 ) = ‘ lag/l/ellipse/ total ' , 
sourcetype (7) = • lag/l/ellipse/ total ' , 
sourcetype < 8 ) = ' lag/1/ ellipse/ total ' , 
sourcetype (9)=' lag/l/ellipse/ total ' , 
sourcetype ( 10) = * lag/ 1/ ell ipse/ total * , 
sourcetype (11) = ' lag/1/ ellipse/ total ' , 
sourcetype ( 12 ) = ' lag/1/ ellipse/ total ' , 
sourcetype (13) = ' lag/l/ellipse/ total ' , 
sourcetype (14 ) = ’ lag/l/ellipse/ total ' , 
sourcetype ( 15 ) = *lag/l/ellipse/total' , 
sourcetype (16) = * lag/1/ ellipse/ total ' , 
sourcetype ( 17) = * lag/l/ellipse/ total ' , 
sourcetype (18) = ' lag/l/ellipse/ total * , 
sourcetype (19) = ’ lag/l/ellipse/ total ' , 
sourcetype (20) = ' lag/l/ellipse/ total ' , 
sourcetype (21) = ' lag/l/ellipse/total ' , 
numparts= 200, 200, 

200, 200, 200, 20( 

200, 200, 200, 20( 

0, 0, 0, 0, 

0, 0, 0, 0, 

0, 0, 0, 0, 

0, 0, 0, 

rotation= 0 . OOOOOOOOOOE+OO , 
nuccg= 400, 

mycg= 400, 

mzcg= 40, 

delxcg= 600.000000 
delycg= 600.000000 
delzcg= 75.00000000 , 

centxcg= -80.56600 , 
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28.00 


centycg= 

ixyslab= 

ixzslab= 

iyzslab= 

&end 


4 , 

100 , 

0 , 



APPENDIX D 

PROWESS HYPACTJN file 
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&HYPACT_IN 

hypref = ' / ul/met/970 117 . 12 / RAMS/ f . a ' , 
f reqopwrt=600 . , iavgout=l/ 

hpartfile=’ ./test/H. ' , 
dtpart=120 , , 

maxpart=46000, npartbl=46000 , 

f reqavg=600 . , 

ihturb=l , 

iadvord=2 , 

nspecies=0 , 

npsources=2, 

irelstrt=163000, 163000, 060000, 060000,060000, 
istrtdays= 0, 0, 0, 0, 0, 

ireldur =000500,000500,002000,002000,002000, 
idurdays= 0, 0, 0, 0, 0, 

isiinend =203000, 
ienddays= 0 , 

nuinparts=3000 , 3000, 400,400, 400, 
emission=l .,1.,1.,1.,1*, 
tricoord= -4. ,42, , -3 . ,43. , -2. ,41. , 

2. ,4. ,6. ,3. ,5. ,7., 

2. ,4. ,6. ,3. ,5, ,7. , 
srcy=28 .448,28,448, 
srcx= -80.566, - 80.566, 
srcz=600 . , 1800 , , 
xsize=600 . ,600. ,3000. ,3000. , 
ysize=600 . ,600. , 6000. , 6000. , 
zsize= 4 00 . , 400 , , 100 . , 100 . , 
source type ( 1 ) = ' lag/ 1/ ellipse/ total * , 

sourcetype (2) = ' lag/ 1/ellipse/ total' , 

source type (3) = ' lag/1/ ellipse/ total ’ , 
sourcetype ( 4 ) = ' lag/ 2/ rectangle/ total , 

sourcetype (5) = *both/3/ellipse/scaled' , 

rotation=0 . , 0 . , 30 . , 30 . , 45 , , 45 . , 
ihfall=0, 

szpwr=l .,l.,l.,l-,l.r 
szmin=l -,1.,1.,1.,1*» 
szmax=l .,1.,1.,1./1-/ 
mxcg=400 , 
mycg=400 , 
mzcg=4 0 , 
delxcg=600 . , 
delycg=600 . , 
delzcg=75 . , 
centycg=28 . 00 , 
centxcg= - 80 . 566 , 
ixyslab=4 , 
ixzslab=100 , 
iyzslab=0 , 

&end 



NOTICE 


Mention of a copyrighted, trademarked or proprietary product, service, or document does not constitute 
endorsement thereof by the author, ENSCO, Inc., the AMU, the National Aeronautics and Space 
Admimstration, or the United States Government. Any such mention is solely to inform the reader of the 
resources used to conduct the work reported herein. 
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